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survey,  conducting  the  survey  and  reporting  on  the  survey  were  included  in  the 
study, 

Tl$e\site  modeling  program  models  the  manufacturing  processes  which  contri- 
bute pollutants  to  the  system,  models  the  sewer  system,  and  models  the  treat- 
ment system\including  acid  or  caustic  neutralization,  settling  ponds,  and 
domestic  treatment.  The  inputs  to  the  model  are  the  production  levels  of  the 
manufacturing  processes  and  the  outputs  are  the  predicted  pollutant  measurement 
values  at  each  possible  measure  point  in  the  system. 

The  resource  matching  program  accepts  data  defining  proposed  measurements 
and  matches  these  against  the  available  time,  manpower,  and  equipment.  The 
output  lists  the  pollutant  to  be  measured  at  each  measure  point,  the  total 
commitment  of  time  for  each  analyst  and  for  each  piece  of  equipment.  Note  is 
made  of  any  overcommitment  of  manpower  or  equipment. 

The  model  refinement  or  updating  program  accepts  measurements  taken  durin 
a preliminary  survey  or  during  a regular  survey  and  computes  suggested  new 
parameters  for  the  process  models. 

The  indicator  model  program  evaluates  the  performance  of  sanitary  treat- 
ment facilities. 

The  program  uses  design  data,  data  from  the  operating  log  and/or  data 
generated  during  the  survey  and  computes  key  operational  characteristics.  Com- 
paring these  with  desirable  values  as  cited  in  design  books  and  manuals  will 
give  the  survey  planner  insight  into  the  operation  of  the  system  and  suggest 
the  need  for  more  survey  measurements  or  the  need  for  changes  in  operation. 

A system  was  developed  for  automatic  instrumentation  of  pH,  conductivity, 
and  other  parameters  which  use  strip  charp  recordings.  Interface  hardware  was 
selected  and  purchased  and  interface  software  was  developed  for  direct  connec- 
tion to  a digital  computer. 

A data  handling  system  was  developed  for  use  during  and  after  the  survey. 
A PDP8-0S/8  and  peripheral  equipment  was  purchased.  Software  was  developed  to 
perform  data  handling  functions  and  to  direct  the  user  in  application  of  the 
software.  The  program  accepts  raw  data  from  the  analytical  chemist  and  per- 
forms data  conversions,  transcriptions,  and  data  logging  functions.  Output 
is  available  in  several  forms  as  may  be  needed  for  various  reports  during  and 
at  the  end  of  the  survey. 

Recommendations  are:  the  survey  planner  should  obtain  sufficient  data  in 

a preliminary  survey  to  model  and  analyze  the  site;  measurements  should  'e  aut 
mated  to  the  maximum  extent  possible;  data  handling  should  be  delegated  to  tne 
computer  when  the  operations  are  well  defined  and  repetitive.  The  programs, 
software  and  hardware  included  here  will  assist  the  survey  planner  in  followin 
these  recommendations  and  design  a more  effective  survey. 
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ABSTRACT 

This  is  the  final  report  of  a three  year  project  titled,  "A  Systems  Analysis 
of  Water  Quality  Survey  Design." 

In  this  project  a study  was  made  of  water  quality  surveys  conducted  by  the 
United  States  Army  Envl ronmental  Hygiene  Agency  (AEHA) . Mainly  data  and  reports 
from  studies  of  Army  Ammunition  Plants  (AAP)  were  used. 

The  focus  of  this  project  was  the  development  of  computer  aided  procedures 
which  would  assure  efficient  use  of  manpower  and  equipment  and  assure  that  the 
measurements  taken  give  a reasonable  representation  of  the  system.  Planning  the 
survey,  conducting  the  survey  and  reporting  on  the  survey  were  included  in  the 
study. 

The  site  modeling  program  models  the  manufacturing  processes  which  contri- 
bute pollutants  to  the  system,  models  the  sewer  system,  and  models  the  treatment 

« 

system  including  acid  or  caustic  neutralization,  settling  ponds,  and  domestic 
treatment.  The  inputs  to  the  model  are  the  production  levels  of  the  manufactur- 
ing processes  and  the  outputs  are  the  predicted  pollutant  measurement  values  at 
each  possible  measure  point  in  the  system. 

The  resource  matching  program  accepts  data  defining  proposed  measurements 
and  matches  these  against  the  available  time,  manpower,  and  equipment.  The  out- 
put lists  the  pollutant  to  be  measured  at  each  measure  point,  the  total  commit- 
ment of  time  for  each  analyst  and  for  each  piece  of  equipment.  Note  is  made  of 
any  overcommitment  of  manpower  or  equipment. 

The  model  refinement  or  updating  program  accepts  measurements  taken  during 
a preliminary  survey  or  during  a regular  survey  and  computes  suggested  new  para- 
meters for  the  process  models. 

The  indicator  model  program  evaluates  the  performance  of  sanitary  treatment 
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The  program  uses  design  data,  data  from  the  operating  log  and/or  data  generated 
during  the  survey  and  computes  key  operational  characteristics.  Comparing 
these  with  desirable  values  as  cited  in  design  books  and  manuals  will  give  the 
survey  planner  insight  into  the  operation  of  the  system  and  suggest  the  need 
for  more  survey  measurements  or  the  need  for  changes  in  operation. 

A system  was  developed  for  automatic  instrumentation  of  pH,  conductivity, 
and  other  parameters  which  use  strip  chart  recordings.  Interface  hardware  was 
selected  and  purchased  and  interface  software  was  developed  for  direct  connection 
to  a digital  computer. 

A data  handling  system  was  developed  for  use  during  and  after  the  survey. 

A PDP8-OS/8  and  peripheral  equipment  was  purchased.  Software  was  developed  to 
perform  data  handling  functions  and  to  direct  the  user  in  application  of  the 
software.  The  program  accepts  raw  data  from  the  analytical  chemist  and  performs 
data  conversions,  transcriptions,  and  data  logging  functions.  Output  is  avail- 
able in  several  forms  as  may  be  needed  for  various  reports  during  and  at  the  end 
of  the  survey. 

Recommendations  are:  the  survey  planner  should  obtain  sufficient  data  in 

a preliminary  survey  to  model  and  analyze  the  site;  measurements  should  be  auto- 
mated to  the  maximum  extent  possible;  data  handling  should  be  delegated  to  the 
computer  when  the  operations  are  well  defined  and  repetitive.  The  programs, 
software  and  hardware  included  here  will  assist  the  survey  planner  in  following 
these  recommendations  and  design  a more  effective  survey. 
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I.  Introduction 

This  research  addresses  the  problem  of  improving  the  design  of  short-term 
(15  days)  water  quality  surveys.  These  surveys  are  primarily  concerned  with 
monitoring  for  industrial  wastes  generated  at  Army-owned  and,  sometimes,  con- 
t ractor-operated  facilities.  Most  of  the  data  used  was  from  munitions  facil- 
ities. 

In  general  terms  a water  quality  survey  involves  a survey  team  going  to 
the  site  (e.g.  an  ammunition  plant)  and  over  a period  of  several  days  collecting 
water  samples.  These  samples  are  collected  at  predetermined  time  intervals  at 
predetermi ned  points  in  the  site  sewer.  The  frequency  of  sampling  and  points 
for  sampling  are  part  of  a survey  plan  or  design  determined  after  a preliminary 
site  visit.  The  samples  are  analyzed  to  determine  the  level  of  pertinent 
pollutants  in  the  samples.  Some  of  the  analyses  are  done  on  site  in  portable 
laboratory  vans  and  some  of  the  analyses  are  performed  in  a central  laboratory 
facility.  This  latter  requires  sample  preparation  and  shipping.  A survey 
report  is  prepared.  The  report  clearly  defines  the  sampling  and  analysis 
procedures  used  and  lists  survey  results  with  tables  and  curves.  Thus,  the 
report  identifies  any  time  and  location  variation  in  pollutant  levels  measured. 
Also,  any  potential  pollution  problem  and  its  cause  should  be  detected. 

A good  survey  design  is  one  which  yields  data  which  accurately  describes 
the  system  state  with  a minumum  of  redundant  data.  To  paraphrase  the  previous 
statement,  a good  survey  design  is  one  which  yields  a maximum  of  useful  in- 
formation within  the  limitation  of  the  available  time,  personnel,  and  equipment. 

In  a review  of  reports  of  past  surveys,  it  was  noted  that  some  of  the 
pollutants  measured  at  some  of  the  measure  points  were  at  levels  too  low  to  be 
significant  or  remained  unchanged  throughout  the  survey  period.  If  modeling 
and  preliminary  analysis  indicated  one  of  the  above  cases  to  apply  to  a planned 
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survey,  analysis  for  that  particular  pollutant  on  one  or  two  samples  would  be 
sufficient  to  verify  the  preliminary  analysis  and  much  of  the  chemists  time 
could  be  saved  or  put  to  better  use.  Modeling  and  preliminary  analysis  can 
also  be  used  to  trace  any  pollutant  back  to  its  source  and  this  clearly  defines 
a cause-effect  relationship.  A good  survey  design  should  not  allot  time  for 
measurement  of  i ns  i gn i f icant  parameters  and  the  results  should  point  out  causes 
for  presence  of  significant  pollutants. 

Efficient  use  of  time,  personnel,  and  equipment  is  part  of  a good  survey 
plan.  For  each  sample  collected  several  chemical  analyses  are  performed  to 
determine  the  level  of  pertinent  pollutants.  In  general  each  data  point  involves 
a chemical  analysis  and  not  merely  a meter  reading.  Tests  for  different 
pollutants  may  require  different  equipment,  different  amounts  of  time,  and 
different  abilities  for  the  analyst.  A good  survey  design  should  include  a 
correlation  between  the  number  of  samples,  the  number  of  analyses  on  each  sample, 
and  the  time,  personnel  and  equipment  available  for  the  survey. 

Good  data  management  should  be  a corollary  to  a good  survey  plan  or  design. 
As  a general  rule,  any  data  manipulation  that  is  well  defined  and  highly 
repetitive  should  be  delegated  to  the  digital  computer.  When  hardware  is 
available,  analysis  should  be  automated  and  interfaced  with  the  digital  computer. 
When  automation  is  not  practical,  raw  analysis  data  should  be  fed  to  the  computer 
by  the  analytical  chemist  and  the  computer  should  be  programmed  to  perform  con- 
versions, transcriptions  or  logging.  Preparation  of  data  for  reports,  plotting 
curves,  and  setting  up  tables  are  valid  computer  operations.  Reports  are  needed 
during  and  at  the  end  of  a survey.  The  computer  should  be  programmed  to  generate 
these  reports  or  appropriate  output  for  these  reports. 

The  design  of  water  quality  surveys  was  invertlgated  by  an  interdisciplinary 
team  using  systems  analysis  techniques.  The  team  Included  faculty  and  graduate 
student  researchers  with  expertise  in  environmental  systems,  environmental 
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chemistry,  computer  programming,  computer  hardware  and  systems  analysis.  As 
the  research  progressed,  six  objectives  evolved.  Work  on  these  objectives 
i ncl uded : 

1.  A method  was  developed  for  computer  modeling  of  the  site  to  be  surveyed. 
The  various  manufacturing  processes  of  the  site  form  the  pollutant  sources.  Each 
process  model  has  as  its  input  the  production  level  and  as  its  output  the 
effluent  it  discharges  into  the  sewer  systems.  The  effluent  is  defined  in  terms 
of  flow  rate  and  a parameter*  measure  number  for  each  pertinent  pollutant.  The 
sewer  system  into  which  the  processes  discharge  their  effluent  is  modeled  such 
that  for  any  branch  of  the  sewer  system  the  contribution  of  each  source  is 
determined  and  hence  to  total  flow  and  each  pollutant  parameter  for  that  branch. 

The  survey  planner  can  hus  specify  the  production  level  of  each  process  on  the 
site  and  determine  the  flow  and  pollutant  parameters  for  each  branch  of  the  sewer 
system.  The  site  model  also  allows  for  caustic  or  acid  neutralization,  settling 
ponds  and  domestic  treatment  within  the  sewer  system.  By  specifying  the  production 
levels  that  may  occur  for  each  site  process,  the  survey  planner  is  better  informed 
in  deciding  what  should  be  measured,  where  it  should  be  measured  and  how  often 

it  should  be  measured. 

2.  A procedure  was  defined  for  modeling  the  measurement  resources.  Also 
a program  was  developed  which  utilizes  the  site  model  and  the  measurement  re- 
sources model  and  determines  the  commitment  of  resources  required  for  a measure- 
ment plan  proposed  by  the  survey  planner.  The  resource  model  lists  personnel, 
equipment,  and  time  required  to  measure  each  parameter  for  each  method  available 


The  term  parameter  measure  number  or  just  parameter  is  used  here  since  some 
pollutants  are  measured  in  terms  of  concentrations  such  as  milligrams  per 
liter  while  others  such  as  pH  and  turbidity  each  have  their  own  units. 
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The  output  of  the  program  lists  the  method  of  measurement  of  each  parameter  to 
be  measured  at  each  point  where  it  is  to  be  measured.  Equipment  used  is  also 
listed.  Total  time  demand  on  personnel  and  equipment  is  given  and  this  is 
compared  with  time  available  for  personnel  and  equipment.  Non-feasible  condi- 
tions are  noted.  By  iteration  through  several  plans  the  survey  planner  can 
arrive  at  a "best"  survey,  a plan  which  yields  a maximum  of  useful  information  for 
the  available  time,  personnel,  and  equipment. 

In  choosing  a survey  plan  the  planner  can  determine  which  pollutants  will 
be  measured  at  which  points  in  several  ways. 

a.  The  program  logic  selects  the  minimum  number  of  measure  points  required 
to  detect  all  of  the  pollutant  arriving  at  the  outfall.  Low  concentrations  may 
preclude  measurement  at  the  outfall. 

b.  The  program  can  be  made  to  select  as  a measure  point  any  source  which 
contributes  more  than  X%  to  the  output  mass  of  a particular  parameter. 

c.  Any  parameter  at  any  point  can  be  manually  flagged  as  a measure  point. 

3.  A program  was  developed  which  will  use  the  site  model  data  and  survey 

results  and  prepare  data  for  verifying  or  updating  the  site  model.  Experience 
in  the  three  years  of  the  study  showed  that  good  dependable  process  data  was 
elusive.  Further,  some  variance  was  noted  in  the  same  processes  at  different 
sites  and  in  different  seasons.  Over  a period  of  time,  use  of  this  program  will 

result  in  more  dependable  process  models.  With  confidence  in  the  process  models, 
the  site  model  can  be  used  to  extrapolate  and  determine  pollutant  conditions  for 
operating  levels  that  exist  at  times  when  measurements  are  not  taken. 

A.  A program  was  developed  for  sanitary  treatment  facilities  (called  the  in- 
dicator model)  which  will  calculate  certain  parameters,  generally  termed  process 
variables  or  operating  characteristics,  that  Indicate  the  operational  state  of  the 
system.  These  calculated  values  can  then  be  compared  with  the  desirable  values  as 
cited  in  design  books  or  manuals. 
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A comparison  will  give  the  engineer  insight  into  the  operation  of  the  system  such 
that  problem  areas  and  possible  causes  of  these  problems  are  more  readily  apparent. 

5.  The  feasibility  of  using  automated  measurement  equipment  was  investigated. 
Direct  measurement  probes  and  interface  (hardware  and  software)  with  the  computer 
were  investigated.  Equipment  was  purchased  and  software  was  developed  for  pH 

and  conductivity  measurements.  A digitizer  was  purchased  and  interfaced  with 
the  computer  for  digitization  of  strip  chart  recordings.  There  is  limited  auto- 
matic measuring  euipment  available  for  water  quality  analysis. 

6.  A data  handling  system  was  developed.  A Digital  Equipment  Corporation 
PDP-8  computer  system  with  OS/8  operating  system  software  and  peripheral  equations 
was  purchased.  Software  was  developed  to  perform  the  data  handling  function  and 
direct  the  user  in  application  of  the  software.  The  program  accepts  raw  data 
from  the  analytical  chemist  and  performs  data  conversions,  transcriptions,  and 
data  logging  functions.  Maximum,  minimum,  mean,  variance,  and  standard  deviations 
are  computed.  Ordering,  classifying  and  storing  operations  are  performed.  At  the 
request  of  the  user,  output  is  available  in  several  forms  as  may  be  appropriate 
for  the  various  reports  needed  during  the  survey  and  at  the  end  of  the  survey. 

These  forms  include  tables  of  data  and  plotted  curves. 

The  philosophy  and  structure  of  the  models,  programs,  and  software  is  described 
in  subsequent  sections  in  the  body  of  this  final  report.  Appendices  to  this  report 
include  users'  manuals,  programmers'  manuals,  program  listings,  and  software  tapes. 
The  manuals  give  the  details  of  procedures,  data  required,  and  options  available 
in  using  programs,  models  and  software.  The  appendices  include: 

A.  The  Survey  Planning  Program 

a.  program  listing  with  example  results  (Appendix  l) 

b.  user's  manual  (Appendix  II) 

c.  programmer's  manual  (Appendix  III) 
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The  card  decks  are  not  included  in  the  appendix  but  were  presented  to  AEHA. 

| i 

The  decks  include: 

Deck  I.  Contains  MAIN  and  subroutines  which  will  permit  modeling 
of  plant  processes  and  sewer  systems.  The  output  includes  the  pol- 
lutant level  at  each  possible  measure  point  in  the  system  for  each 
pollutant  in  the  study.  The  program  allows  consideration  of  acid 
or  caustic  neutralization,  settling  ponds,  and  domestic  sewage  treat- 
ment. Process  production  levels  can  be  varied  for  senstivity  ana- 
lysis. 

Deck  2.  Contains  all  of  Deck  1 but  also  includes  subroutines  which 
will  permit  modeling  of  measurement  resources  and  subroutines  which 
determine  the  commitment  of  resources  for  a measurement  plan  pro- 
posed by  the  survey  planner.  By  Iteration  the  planner  can  arrive 
at  a best  plan  for  available  resources. 

Deck  3-  Contains  all  of  Deck  1 but  also  includes  a subroutine  CORRCT 
which  compares  site  model  data  and  survey  results  and  prepares  data 
for  verifying  the  model  or  updating  the  process  models. 

B.  Domestic  Treatment  Evaluation  (INDICATOR  MODEL) 

a.  program  listing  (Appendix  IV) 

b.  user's  manual  (Appendix  IV) 

c.  programmers'  manual  (Appendix  IV) 

This  program  computes  key  operational  characteristics  which  can  be  compared 
to  desirable  values  for  domestic  treatment  systems. 

C.  Automated  Instrumentation 

a.  user's  manual  (Appendix  V) 

b.  software  to  interface  with  computer  (delivered  to  AEHA) 

c.  programmers'  manual  (Appendix  VI) 
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D . Data  Hand) i ng 

a.  user's  manual  (Appendix  VI  I ) 

b.  software  (delivered  to  AEHA) 

c.  programmers'  manual  (Appendix  VII I ) 

The  software  developed  is  essentially  interactive,  i.e.  it  prints  messages 
which  tell  the  user  what  is  required  next  with  minimum  use  of  the  manual.  This 
system  consists  of  a PDP8-0S/8  and  software  developed  for  accepting  raw  measure- 
ment data,  for  preparing  data  for  reports,  and  for  generating  reports. 

Recommendations  as  a result  of  these  studies  are  that  the  survey  planner 
should  obtain  sufficient  data  in  a preliminary  survey  to  model  and  analyze  the 
site,  that  measurements  should  be  automated  to  the  maximum  extent  possible  with 
commercially  available  equipment,  and  that  data  handling  should  be  delegated  to 
the  computer  at  all  points  in  the  survey  where  the  operations  on  data  are  well 
defined  and  repetitive.  If  the  survey  planner  has  a more  complete  picture  of  the 
site  through  modeling,  redundant  or  meaningless  measurements  could  be  eliminated 
and  more  useful  data  could  be  obtained  in  the  limited  time  available.  In 
automating  measurements  and  in  delegating  repetitive  tasks  to  the  computer, 
efficiency  is  improved.  Further,  requiring  creative  personnel  to  perform 
repetitive  tasks  can  be  counter-productive. 
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II.  SURVEY  PLANNING  PROGRAMS 

A.  Site  Modeling  and  Simulation 

The  purpose  of  site  modeling  and  simulation  In  survey  planning  is  to 
give  the  survey  planner  a detailed  and  accurate  picture  of  the  plant  site 
so  that  he  will  be  better  informed  when  deciding  what  is  to  be  measured  and 
where  the  sample  points  should  be  chosen.  The  preliminary  survey  or  site 
visit  should  give  sufficient  information  to  model  the  plant  and  should  give 
the  survey  planner  information  on  levels  of  production  of  various  parts  of 
the  plant  expected  during  the  survey  period.  Modeling  and  simulation  should 
give  a reasonable  estimate  of  the  quantities  to  be  measured. 


SITE  MODEL  COMPONENTS 


Figure  1 


1 1 

The  components  of  the  plant  site  model  are  shown  in  Figure  1.  The 
components  are:  the  process  model  which  represents  the  plant  operation; 

the  treatment  model  which  represents  the  treatments  that  take  place  within 
the  plant  site;  the  topology  model  which  models  the  plant  site  sewer  system 
and  with  the  process  and  treatment  data  determines  the  state  of  the  effluent 
in  each  branch  of  the  sewer  system  (this  is  output);  the  survey  officer  who 
observes  the  system  states  and  makes  appropriate  changes  in  the  plant  process 
levels  to  represent  expected  plant  production  schedules. 

The  following  sections  describe  these  component  models  in  detail. 

A . 1.0  Process  Model: 

This  model  is  designed  to  define  a manufactur ing  plant  in  terms  of  the 
pollutants  produced  by  the  manufacturing  operations  as  a function  of 
production  levels  of  the  various  plant  components.  In  order  that  the  model 
be  applicable  to  more  than  one  plant,  plants  are  subdivided  on  the  basis  of 
objectives  and  then  further  subdivided  into  processes  required  by  these 
objectives.  A given  plant  may  have  one  or  several  objectives.  Examples  of 
objectives  are:  continuous  TNT;  TNT  (batch);  paint  stripping  of  artillery 

shells.  The  components  representing  the  plant  then  are  these  processes 
required  by  a plant  objective.  The  processes  are  manufacturing  operations 
which  are  definable  in  terms  of  production  level  and  resultant  pollutants. 
These  processes  may  be  components  for  one  or  several  plant  objectives  and 
may  occur  in  one  or  several  plants.  A plant  process  model  is  obtained  by 
assembling  the  appropriate  component  process  models.  The  input  to  the  plant 
model  is  the  production  level  of  each  of  the  processes  making  up  the  plant 
model.  The  output  is  an  ordered  listing  of  the  types  and  levels  of  pollutants 
generated  by  producing  products  at  that  level.  Also  listed  are  the  points 
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in  the  effluent  drain  system  at  which  the  pollutants  enter.  Any  time 
variation  in  pollutant  flow  requires  a new  listing,  a listing  for  each 
state.  The  modeling  program  PROCES  is  written  so  that  the  user  inputs 
to  the  program  only:  the  number  of  processes  at  the  site;  the  process 

codes  for  all  processes  at  the  site  (partial  list  given  in  Table  11, 
complete  list  in  user's  manual);  the  production  level  for  each  process. 
"Process"  models  have  been  constructed  for  eight  plant  "objectives". 

These  objectives  are:  Continuous  TNT;  Nitrocellulose  (batch);  Continuous 

Nitroglycerine;  Single-Base  Power;  "Detrex"  Preparation  of  Forgings; 

Paint  Stripping  of  Artillery  Shells;  Motor  Pool  Wash  Racks;  Troop 
Installation  (sanitary).  All  of  the  process  models  are  indentified  in  a 
list  in  Table  l of  the  user's  manual  in  the  Appendix.  The  process  models 
are  given  in  the  Appendix  in  the  Program  Listing. 

The  pollution  potential  of  an  effluent  is  measured  in  a number  of 
ways.  One  type  of  measurement  is  to  determine  the  concentration  of  a 
particular  element  or  compound.  An  example  of  this  type  of  measurement 
would  be  the  measure  of  sulfates  in  milligrams  per  liter.  Another  type  of 
measurement  measures  a character i st ic  of  the  system  not  necessarily  related 
to  only  one  chemical  compound  of  element.  Examples  of  this  type  of  measure- 
ment are  conductivity,  color,  and  pH.  To  speak  in  general  terms  of  the 
numbers  that  define  the  level  of  pollution  potenital,  whether  they  are  of  one 
or  the  other  of  the  above  types,  the  term,  "pollution  parameter"  will  be  used. 

Figure  II  is  an  example  of  a computer  program  subroutine  model  of  a 
process.  All  of  the  other  processes  are  given  in  the  appendix  in  the  Program 
Listing.  Statement  0001  gives  the  subroutine  name.  Statement  0002  is  the 
dimension  statement.  Statement  0003  is: 
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] . pH 

2.  COD 

3.  CONDUCTIVITY 
k.  TOC 

5.  AMMONIA  (NH3)  - N 

6.  DISSOLVED  SOLIDS 

7.  NITRITE/NITRATE 

8.  TOTAL  N (KJELDAHL) 

9.  TOTAL  SOLIDS 

10.  SUSPENDED  SOLIDS 

11.  TOTAL  HANDNESS 

12.  SULFATES 

13.  MBAS 

14.  SULFITES 

15.  CHLORIDES 

16.  VOL.  SUSP.  SOLIDS 
17-  SODIUM 

18.  TOTAL  PHOSPHATES 
19-  TNT 

20.  CALCIUM 

21.  COLOR 

22.  ALKALINITY 

23.  TURBIDITY 
2k.  ACIDITY 
25-  BOD 

MASTER  LIST  OF  PARAMETERS 
Table  1 
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EXAMPLE  PROCESS  MODEL 
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FLOW ( JM) =0. 01 *CAP ( JM) *0. 10 

This  statement  computes  the  effluent  flow  from  process  numbered  (JM)  in 
million  gallons  per  day  (MGD)  for  a given  operating  level,  CAP(JM)  in 
% capacity.  The  0.01  premultiplier  changes  the  $ to  a decimal  fraction 
and  the  0.10  postmultiplier  is  the  MGD  of  the  process  at  1 00%  capacity. 

The  remaining  lines  0004  through  0019  give  the  pollutant  parameters 
necessary  to  define  this  process.  These  parameters  must  be  indentified 
by  a number  assigned  in  a master  list  given  in  Table  1.  For  example, 

SPLIST  (JM,3)  * 4000  gives  the  conductivity  of  the  effluent  for  the 
process  numbered  JM.  Three  is  the  parameter  number  for  conductivity  from 
the  master  list. 

Further,  each  process  is  indentified  by  a mnemonic  from  a master  list. 

A partial  list  is  given  in  Table  11.  A complete  list  is  in  the  User's 
Manual  of  the  Appendix.  The  mnemonic  for  calling  this  subroutine  is  PTNT 
and  the  subroutine  name  is  STNT. 

The  processes  are  numbered  in  the  order  in  which  they  are  read  in  by 
the  process  model.  The  ordering  or  the  numbering  is  arbitrary  except  that 
it  must  be  consistent  with  the  numbering  used  in  the  TOP  subroutine  modeling 
the  sewer  system.  This  data  is  read  in  the  START  subroutine.  It  is  by  this 
numbering  that  the  modeling  disignates  where  a process  discharges  into  the 
sewer  system. 

The  plant  process  modeling  in  review  is  as  follows:  During  the 

preliminary  survey  all  processes  which  make  up  the  plant  operation  are 
identified  in  terms  of  the  proper  mnemonic  from  the  master  list.  The  expected 
operating  capacity  in  I of  full  load  Is  noted  for  each  process.  The  mnemonic 
and  capacity  are  placed  on  a data  card  for  read  In  by  the  PR0CES  subroutine. 
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CORRESPONDING 

SUBROUTINES 

MNEMONIC  - DESCRIPTION 

PAPC  - AOP  Compression  SAPC 

PAOP  - AOP  - Flow  3 SAOP 

PSAC  - Sulfuric  Acid  Concentrator  SSAC 

PNAC  - Nitric  Acid  Concentrator  SNAC 

PBDN  - Blowdown  at  BAAP  SBDN 

PB I X - Ion  Exchange  at  BAAP  SBIX 

PCOW  - Cooling  Water  - Flow  10  SCOW 

PNBP  - Nitration  and  Boiling  Operations  SNBP 

PBPP  - Beater  and  Poaching  Operation  SBPP 

PBEX  - Benzene  Extraction  SBEX 

PNCS  - Nitrocellulose  Into  Ball  Power  SNCS 

PNGS  - Ball  Power  with  Nitroglycerin  SNGS 

PPND  - Pond  at  Head  of  Sewer  SPND 

PAPP  - AOP  - Flow  A SAPP 

PSWG  - Sewer  Influent  - Flow  into  STP  Process  SSWG 

PCWZ  - AAP  - Cooling  Water  - Flow  11  SCWZ 

PPAS  - Purification  Acid  Scrubber  SPAS 

PPCU  - Purification  Clean-Up  SPCU 

PPCW  - Purification  Cooling  Water  SPCW 

PFSW  - Finishing  Scrubber  Water  SFSW 

(Partial  List) 

PROCESSES  AND  MNEMONICS 
(SOURCES) 


Table  II 
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The  process  cards  are  read  in  the  order  defined  later  in  the  discussion  of 
the  TOP  subroutine. 

As  the  subroutine  PROCES  reads  in  the  process  data  cards  it  forms  a 
plant  model  in  terms  of  a SPLIST  matrix.  The  SPUST  matrix  has  as  its  rows 
the  different  processes  present  and  as  its  columns  the  parameter  values, 
the  parameters  in  order  as  defined  by  the  master  list  of  Table  II  and  the 
processes  in  the  order  read  in.  An  abbreviated  SPLIST  matrix  format  for 
five  processes  and  four  pollutants  is  given  in  Figure  IN. 


PH 

COO 

CONDUCTIVITY 

TOC 

Pollutant  Names 
in  Columns 

PPND 

6.0 

20.0 

0.0 

11.0 

PAPC 

6.2 

20.0 

0.0 

11.0 

Pol  1 utant 

PAOP 

6.2 

20.0 

0.0 

11.0 

Parameters 

PAPP 

6.2 

20.0 

0.0 

11.0 

PNAC 

2.5 

25.0 

0.0 

22.0 

Process  Model  Names  in  Rows. 

EXAMPLE  SPLIST  MATRIX 
Figure  I | | 

An  important  feature  of  the  PROCES  subroutine  is  that  it  will  accept 
data  which  will  alter  one  or  more  parameters  for  any  process  present.  If 
the  user  feels  that  at  a particular  plant  site  a different  parameter  value 
would  give  a better  model,  the  parameter  value  can  be  changed.  Further  in 
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a more  extreme  case  all  new  parameter  values  can  be  read  in  for  a given 
process  if  existing  data  and  the  judgement  of  the  user  warrant  the  change. 
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The  procedure  for  making  these  changes  Is  given  in  the  user's  manual. 
Clearly  this  latter  capability  can  be  used  to  Include  a process  in  a 
survey  planning  when  the  process  Is  not  In  the  process  library. 
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A . 2.0  S i te-Topoiogy  Model: 

The  site-topology  model  has  as  its  body  a codified  listing  of  the  details 
of  the  physical  layout  of  the  plant  site  sewer  system  and  a program  which 
computes  the  flow  and  pollution  parameter  values  for  each  branch  of  the  system. 

The  details  of  the  physical  layout  of  the  plant  site  sewer  system 
includes: 

1.  I ndent if icat ion  of  all  processes  (sources  of  pollution)  and  the 
location  at  which  they  discharge  into  the  sewer  system. 

2.  I ndent i f icat ion  of  all  manholes,  weirs  or  other  possible  measure 
points . 

3-  I ndent i f icat ion  of  all  processes  effluents  which  flow  through  each 
of  the  possible  measure  points. 

When  the  system  is  large  it  is  recommended  that  the  recording  and  coding 
of  the  physical  layout  data  be  done  in  three  steps  as  follows: 

1.  During  the  preliminary  survey,  construct  a map  showing  the  processes, 
the  treatments,  the  sewer  systems,  the  manholes  or  other  measure  points,  and 
the  flow  directions.  Figure  IV  shows  a part  of  Badger  AAP  and  illustrates  the 
data  that  should  be  recorded.  The  blocks  show  processes  and  treatments,  the 
lines  with  arrows  show  the  sewer  lines  and  flow  directions,  and  the  small 
circles  show  the  possible  measure  points. 

2.  The  essential  information  from  the  site  map  is  then  included  in  a 
linear  graph  and  the  Items  are  ordered  or  numbered.  In  the  linear  graph  each 
branch  (line)  represents  a possible  measure  point.  The  processes  (pollutant 
sources)  and  treatments  are  shown  also.  The  graph  shows  each  possible  measure 
point  and  the  effluents  which  flow  through  these  points.  The  source  numbers 
are  assigned  arbitrarily  (must  run  consecutively  starting  from  1).  However, 
once  assigned,  the  same  numbering  must  be  used  In  the  process  model,  the  site- 
topology  model  and  other  parts  of  the  program.  The  linear  graph  branches 
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Ni trocel lulose 
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are  numbered  with  the  following  restrictions: 

a.  The  numbering  must  start  at  1 and  run  consecutively. 

b.  The  numbering  must  be  such  that  in  tracing  from  any  branch  to  the 
outfall,  only  higher  numbered  branches  will  be  encountered. 

c.  Thus,  the  outfall  branch  must  have  the  highest  number.  Figure  V 
illustrates  with  the  linear  graph  for  Figure  IV. 

3.  A convenient  way  of  tabulating  the  topological  data  for  entering 

into  a computer  program  is  to  put  the  information  into  a matrix.  A matrix  is 

constructed  with  a column,  J,  for  each  branch  of  the  system  and  a row,  i,  for 

♦ h 

each  process  of  the  system.  The  entries  for  the  i — row  are  obtained  by 
starting  at  the  i —process  on  the  linear  graph  and  tracing  a path  (with 
the  flow)  through  the  branches  to  the  outfall.  If  the  path  does  not  include 
the  j — branch  the  i j entry  is  zero.  If  the  path  does  include  the  j — 
branch,  an  appropriate  number  is  entered  in  the  i j position.  That  appropriate 
number  depends  upon  whether  or  not  the  flow  undergoes  treatment  in  that  branch 
and  specifically  which  treatment.  If  the  flow  of  source  i passes  through 
branch  j and  no  treatment  occurs  at  branch  j,  then  the  i j entry  is  1.  If  the 
flow  of  source  i passes  through  branch  j and  the  composite  flow  of  branch  j is 
treated  at  branch  j,  then  the  i j entry  is  2,  3 or  A depending  on  the  type  of 
treatment  as  follows: 

Acid  or  caustic  neutralization  2 

settl ing  pond  3 

domestic  sewage  treatment  A 

The  matrix  for  the  linear  graph  of  Figure  V is  shown  In  Figure  VI. 

If  the  sewer  system  is  in  separate  parts  and  hence  has  multiple  outfalls, 
the  outfall  elements  of  the  separate  parts  are  connected  together  and  a 
fictitious  new  outfall  element  is  included  In  the  linear  graph.  In  the  Input 
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POND  (1) 
PARC  (2) 
PAOP  (3) 
PSAC  (4) 
PNAC  (5) 
PAPP  (6) 
PBDN  (7) 
PBIX  (8) 
PCON  (9) 
PCWZ  (10) 


TOPOLOGY  MATRIX  OF  FIGURE  V 
Figure  VI 

data,  a data  item  designated  NTOP  gives  the  program  the  information  necessary 
to  handle  this  situation.  Figure  Vll  illustrates  the  addition  of  the 
f ict i t ious  outfa 1 1 . 


21  - outfal 1 of  Part  I 

22  - outfal 1 of  Part  2 

23  • outfa 11  of  Part  3 

24  - fictitious  outfall 


Figure  VI I 


The  operation  of  the  program  which  computes  the  flow  and  pollutant 
parameter  values  is  illustrated  in  the  following  example. 

The  output  data  of  the  process  models  are  the  inputs  to  the  topology 
model.  The  format  used  for  the  input  data  is  in  vector  form  illustrated  by 
1=1.  (V,,  P.,  p2  — ) 

Where  I.  is  input  at  point  i in  system 
V.  is  volume  flow  rate  of  input 

Pj  is  measurement  level  of  pollutant  parameter  for  the  jth 
parameter . 

The  contribution  of  each  input  to  each  branch  (measure  point)  can  be 
determined  by  well-known  network  topology  theory.  A simple  example  using 
the  system  of  Figure  VIII  will  illustrate: 
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Let  the  given  tree  graph  (solid  lines)  show  the  branches  with  vertices 
lettered  and  the  elements  numbered.  Consider  the  case  where  the  effluent 
waste  lA,  l0>  and  1^  are  entering  the  system  at  A,  B,  and  E,  respectively. 
The  dashed  lines  represent  sources. 

Trace  circuits  In  the  diagram  of  Figure  Vlll  such  that  one  and  only  one 
source  element  Is  Included  In  each  circuit.  Continue  until  all  sources  have 
been  included.  The  circuit  Information  is  tabulated  in  matrix  form  where 
the  columns  correspond  to  the  system  elements  and  the  rows  correspond  to 
the  sources. 
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Also,  let  the  source  matrix  of  elements  6,  7,  and  8 be 


Then  one  can  determine  the  contribution  of  each  source  at  each  section 


of  the  system  by: 


The  (+)  in  the  expressions  1^  + lg  and  1^  + lg  + 1^  is  interpreted  to 
mean  "the  combination  of."  This  combination  procedure  used  here  will  be  a 
direct  sum  for  the  flow  rates  and  a flow  weighted  combination  for  pollutant 
parameter  levels.  The  following  example  illustrates  this: 

Consider  lA  (VA>  P(  ) and  l0  (V0,  P(  ) and  lc  - + lg 


where  VA  - 1,000  gal/hr 

V0  ■>  500  gal/hr 
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P.  * 50  mg/liter  of  sulfate 

'a 


P,  * 75  mg/liter  of  sulfate 

B 

mg  sulfate/hr  of  A » 50  x 1,000  x K 


mg  sulfate/hr  of  B ■ 75  x 500  x K 
K converts  gal  to  liters 


VA  + tfB 


1,000  + 500 


58.4  mg/1 i ters 


thus  lc  (Vc,  P > 


lA  (1000,  50)  + lB  (500,  75) 
lc  (1500,  58.4) 


It  should  be  emphasized  that  the  so-called  "flow  weighting"  defined  in 
the  above  example  is  a first  order  modeling.  It  should  give  reasonably 
accurate  representation  of  combinations  when  considering  such  mass  quantities 
as  milligrams  per  liter  of  sulfates.  However,  for  pH  and  when  treatments  take 
place  more  detailed  calculations  are  performed. 


The  pH  resulting  from  the  mixing  of  the  flows  from  several  sources  is 
calculated  assuming  that  the  alkalinity  and  the  capacity  factor  (CT)  of  the 
combined  flow  are  conservative,  i.e.  can  be  calculated  by  the  above  flow 
weighting  procedures  and  that  the  combined  pH  can  be  calculated  from  the 
alkalinity  and  capacity  factor  of  the  combined  flow.  The  equations  and 
procedures  used  follows: 

Given  the  alkalinity  (ALK)  and  pH  of  each  of  the  incoming  flows,  [H+] 
is  computed  for  each  flow  by 

[H+]  - 10  pM,  the  H-Ion  concentration. 


) 
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The  CT-j.  and  the  first  approximation  of  [H+]  are  used  to  calculate 
an  ALK.  By  Iteration  with  Increments  or  decrements  on  [H+]  the 
calculation  Is  continued  until  an  [H+]  Is  found  which  gives  AUCp.  This 
[H+]  Is  the  value  taken  to  determine  the  combined  flow  pH. 

The  TOP  subroutine  computes  the  flow  and  pollution  parameter  values 
for  each  branch  of  the  system.  It  uses  the  above  illustrated  flow  weighting 
and  pH  calculations  for  each  branch  where  no  treatment  Is  performed.  In 
branches  where  treatment  takes  place,  TOP  calls  In  the  appropriate  sub- 
routine which  modifies  the  branch  parameters  properly  for  the  particular 
treatment  which  takes  place.  The  treatments  and  subroutines  available  in 
this  study  are:  NEUTRA,  acid  or  caustic  neutralization;  SETTLE,  clarifier 

or  settling  pond;  DOME,  domestic  or  sanitary  treatment. 

The  NEUTRA  model  allows  consideration  of  acid  neutralization  using 
Ca(0H)2,  CaCO^,  Na2  CO^,  or  NaOH.  Also,  caustic  neutra 1 izat Ion  is  permitted 
using  HjSO^,  HCI,  HNO^,  or  COj.  The  basic  relations  and  equations  that 
determined  the  model  are  given  in  the  following  paragraphs. 
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A.  2.1  NEUTRA  MODEL 

( 1 )  Analytical  Determinations 

(A)  Ak  la  1 i n i ty 

(1)  [CO.^  - Aik]  = Phenolphthalein  Alkalinity  (eq/£) 

(2)  [Aik]  = Methyl  Orange  Alkalinity  (eq/J.) 

(B)  Acidity 

(1)  [H+  - ACY]  = Methyl  Orange  Acidity  (eq/8.) 

(2)  [CO^  ~ ACY]  = Phenolphthalein  Acidity  (eq/2.) 

(3)  [ACY]  = Titration  to  pure  CO^  = (eq /i) 

(Note:  [ACY]  is  not  routinely  determined 

analytically,  but  is  obtained  by 
calculation. ) 

(11)  Basic  Capacity  Equations 

(A)  [ALK]  = CT  (a,  + 2«2)  + (10'1V[H+]  - [H+] 

(B)  [ALK]  = CT  (2aQ  + a,)  + [H+]  - ( 1 o‘ ' V [H+] ) 

where , 


[H+  1 « do)'pH 

- J 
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where, 

K = i* . i*i*6  X I0"7  g 25°C 
K2  - 4.688  X 10'"  g 25°C 
KfK2  = 2.0843  X lo'17  g 25°C 

III.  Alkalinity/Acidity  Relationships 

A.  [ALK]  + [H-ACY]  - 0 

B.  [ACY]  + [OH  - ALK]  - 0 

C.  [C03  - ALK]  + [C02  - ACY]  - 0 

D.  [ALK]  + [C02  - ACY]  *=  CT 

E.  [ALK]  + [ACY]  + 2 Cy 

F.  [ALK]  - [C03  - ALK]  = Cy 

G.  [C02  - ACY]  - [H-ACY]  = CT 

where,  [OH  - ALK]  = titration  to  pure  C03  * (eq/l) 

IV.  Acid  Neutralization  (H2S0^,  HCI,  HN03)/Using  Ca(0H)2 

(A)  Calculate  [ACY]  in 

(1)  Analytically  determined  [H-ACY]  s [C02  - ACY]  s [pH].n 

(2)  Calculate  (C-j-)jn  using  Equation  lll.G 

(3)  Calculate  [ACY].^  using  Equation  11. B,  (Cy).^,  and  ( pH ) . n 

(B)  Stoichiometric  Equations 

(1)  H2S0^  + Ca(H0)2  ^ CaSO^  + 2H20 

(2)  2HCI  + Ca(0H)2  -*•  CaCl2  + 2H20 

(3)  2HN03  + Ca (OH) 2 * Ca(N03)2  + 2H20 
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Note:  CaCl2  6 Ca(N0^)2  are  very  soluble  while  CaSO^  Is  relatively  Insoluble. 

(C)  Calculate  [Ca (OH) 2]-dd.d 


(1) 

(2) 

(3) 


During  neutralization  with  Ca(0H)2 

*Vout  * ^Vin 

Assuming  neutralization  reaction  is  carried  to 

pH  = 7 ([H+]  - 10  ^),  then  using  Equation 

11. B,  [H+]  * 10  and  (CT)  , calculate  [AC Y ] 

T out  out 

[Ca(0H)2]  added  » ( [AC Y]  in  - [ACY]Qut)/2  (moles/1 i ter) 


(D)  Analytically  Determine  tS04  +] * ; n 

(1)  [SO. ] - 1.04616  X I0'5  * [SO. ] * . 

4J  4J  in 

Moles/liter  mg/£ 


(E)  Calculate  [Cal  in  Moles/Liter 

added  

(l)  [Cal  ad M ‘ Ita(0H)2J»dded  E«"- (l) * *  IV'C3 
moles/liter  moles/liter 


(F)  Calculate  [Ca]  tQta,  jn 

(,)  [Ca]  Total  in  " [Ca  1 + [Ca  ]added 

Analytical  (moles/e.)  Eqn.lV.E.l 
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(H)  If  SOLYPROD  >1 .23  X 10~4  then  have  precipitation  of  CaSO^ 


(1)  [CaSO.l  = Quantity  of  CaSO,  precipitated  (moles/2)  then  ([Ca  ] . 

4 ppt  4 tota I i n 

' I^so„!ppt)«so4  ]in  - !caso4]ppt.  1.32  x .o'* 


(2)  Solve  for  [CaSO,.]  „ (moles/2) 

4 ppt 


+ (1Ca  ■total  it  1SOl’in-  (iCa  'total  ln*tS01  ’J  ^ 

<[Ca  ’tot. I In*150.,  ’in"’-32  * 10 


(3)  [SSI,  - [SSI.  + [CaSO, ] *1.36  * I05 

out  in  4 ppt 

(mg/2.)  (mg/2,) 

(k)  [TDS]  =[TDS] . + [Ca (OH)  ] , *7<t  x 103  - [CaSO,]  t *1.36  x I03 

out  in  Z added  4 ppt 


(mg/2,) 


(5)  [C0ND]QUt  = 1.786  ([TDS]out  -8.8) 


(micromho's/cm) 


(mg/2,) 


(6)  [TS]out  - [TDSlout*  [SS]out  or 

!TS,out*  lTSlln  * 1Ca<0H)2,add,d*7'*  * '°3 
(mg/2,) 


(7)  (pH)out  - 7.0 


(8)  [6LK]out  - 2(CT)out  - [»CY1, 
(eg/2) 


(9)  [Ca]out  - UCa]tou|  |n  - [CaS04]ppt)  *40  x 10- 


TOC  Nitrates 

MBAS  Nitrates 

Su 1 f i tes  Co  lor 

Chlories  Orthophosphate 

Sod i um 

<">  lSVout  - [SOin  - tC*S<Vppt  *96  * '°3 

(12)  Total  Hardness  (TH ) 

[Mg  (mg/X  as  CaCO  )].  = [TH  (mg/X  as  CaCO,)]. 

3 in'  3 i n 

- [Ca  (mg/X  as  CaCO,)]. 

3 i n 

^Mg^out  = [Mg] . 

1 3 i n 

[TH  (mg/X  as  CaCO,)]  . = [Ca  (mg/X  as  CaCO,)] 

3 ^ out  3 out 

+ [Mg  (mg/X  as  CaC0j)].n 


(1)  CSS]  , - CSS]. 

out  i n 

(2)  [TDS]  „ » [TDS] . + [Ca(OH),]  ..  , * 7**  X 103 

out  in  2 added 

(3)  CTS]  ’ = CSS]  + [TDS] 

out  out  out 

or  [TS]  „ - [TS],  + [Ca  (OH)  ] * 71*  X 10 

out  in  2 added 

(mg/Jl) 

(i*)  CC0ND]ou(.  = 1.786  (CTDS]Qut  - 8.8) 

(5)  (pt0out  ■ 7.0 

(6)  [AL«out  - 2(CT)out  - [ACY]out 

(7>  ‘ tCa1in+  ““'added  * W X 11,3 

(8)  lso^-3  - tso4*]|n 

(9)  No  Change:  Same  as  IV.H.1Q 


(10)  Same  formula  as  IV.H.12 


3b 
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V Acid  Neutralization  (H^SO^,  HCI,  HN0^)/Using  CaCO 

(A)  Calculate  [ACY]. 

i n 

( 1 ) Same  as  I V. A. 1 

(2)  Same  as  IV. A. 2 

(3)  Same  as  IV. A. 3 


(2)  Assuming  neutralization  reaction  is  carried  to 
pH  • 7 ([H+]  - 10  ^),  then  using  Equation  II. B, 

[H+]  - 10  ^ and  L ACY 3 , calculate  (C_) 

out  T out 

(3)  [CaC03]added  ■ (CT)out  ' ^ln  <"»»«/>«  ter) 


V 


(D)  Analytically  Determine  [S0^].n* 


( 1 ) Same  as  I V. D. 1 . 


(E)  Ca leu  I ate  [Ca]  added  in  moles/1 i ter 


(1)  [Ca]  = [CaC03laddcd  Eqn.  V.C.3 


>les/£  moles/£. 


(F)  CaJ_c_u]a_te_  [Ca  ]tQta, 
( 1 ) Same  as  I V . F. 1 


(G)  Calculate  Solyprod 

( 1 ) Same  as  I V.G. 1 . 

_ 

(H)  If  Solyprod  > 1.32  x 10  , then  have  precipitation  of  CaSO^ 

(1 ) Same  as  IV. H. 1 

(2)  Same  as  IV. H. 2 

(3)  Same  as  IV.H.3 

(4)  [TDS]out  « [TOS] jn  * tCaC03ladded"  ,o5  - tCaS041ppt*  K36  X 10 

(5)  Same  as  IV. H. 5 

(6)  lTSiout  - [T0S]out  ♦ [SS10U, 


>TS'0U.  ' lTS,ln  * tC»C0J1.d«d*  '°5 
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(7) 

(PH)out  - 7.0 

(8) 

Same  as  IV.H.8 

(9) 

Same  as  IV. H. 9 

(10) 

Same  as  IV.H. 10 

(ID 

Same  as  IV.H. 1 1 

(12) 

Same  as  1 V. H . 1 2 

(1)  If  Solyprod  < 1.32  x 10  ^ 
CaSO^ 

then  have  no  precipitation 

(D 

Same  as  IV.  1 . 1 

(2) 

[TDS]  „ - [TDS] . 
out  in 

* lc*c03>«ld.d*  l<>5 

(3) 

[TS]  _ - [TDS] 
out  out 

or 

+ ^out 

[TS]  - [TS].  + 

out  in 

|CaC03W  '°5 

(*) 

Same  as  IV.  1 . I* 

(5) 

Same  as  IV.  1 .5 

(6) 

Same  as  IV.  1 .6 

(7) 

Same  as  IV. 1 .7 

(8) 

Same  as  1 V . 1 . 8 

(9) 

Same  as  1 V.  1 .9 

(10) 

Same  as  1 V . 1 . 1 0 

Acid  Neutral  1 zet ion  (I^SO^,  HC1 

, HNOj) /Us i ng  Na2C03 

(a)  Calculate  [ACY] . n 


Same  as  IV. A 


J 


p 


1 


(b)  Stoichiometric  Equations 

(1)  H2S04  + Na2C03  - NijSO^  + H2C03 

(2)  2HC1  + N«2C03  -*•  2NaCI  ♦ HjCC^ 

(3)  2HN03  + Na2C03  -*■  2NaN03  + H2C03 


NOTE:  Na„S0,4,  NaCl,  5 NaNO,  are 

very  soluble 

(C) 

Calculate  [Na.CO,] 
2 3 added 

( 1 ) Same  as  V . C . 1 

(2)  Same  as  V.C.2 

<3>  ‘"Wedded  ’ 

<CT>out  • (CT>,n 

(moles/O 

(D) 

‘SS‘out  ' 1551  in 

(E) 

‘Tt,5‘ouc  ‘ 11051  in  * 

‘"Wedded  * ' 

.06  x 105 

(F) 

ITS]  „ - [TDS]  . + 
out  out 

1S5‘ou. 

or  ITS]  „ * [TS] . 

out  in 

+ ‘"W>edd,d  * 

1.06  x I05 

(G) 

[COND]  „ = 1.786  ([TDS]  _ - 8.8) 
out  out 

(H) 

.( PH ) - 7.0 

out 

(1) 

‘AU'‘out  ’ 2(CT>ou.  ' 

■ [ACV‘ou, 

(J) 

‘""‘out  ‘ ‘"“‘‘in  ♦ 

2 ‘‘"Wedded1 

* 23  x 103 

(K) 

NO  CHANGE: 

Same  as  I V.H.  10 
minus  Na 

plus  S0^  6 Ca 


(A)  Calculate  [ACY] 

I n 


Same  as  IV. A 


(B) 


Stoichiometric  Equations 


(1)  H2S04  + 2NaOH  -*•  NajSO^  + 2^0 

(2)  HCI  + NaOH  *-*•  HaCI  ♦ HjO 

(3)  HN03  + NaOH  -*•  NaNOj  + H20 
(C)  Calculate  [NaOH]  ^ ^ 

( 1 ) Same  as  I V. C . 1 

(2)  Same  as  IV.C.2 

(3)  t"a°*ladded  ’ l*"1  - [ACY]  out 

(eq/1  ■ moles/1 ) 

«»  ‘“'out  ' 1551  in 

<t)  lTDS]out  - [TOSj in+  [NaOH]a(Jded  * *0  x 103 

(F)  [TSlout  - (T0S]out  x {SSjQut 

or  [TS]out  • [TS]in+  [»aOH]added  * *0  x I03 

(G)  ICONDJ  = 1.786  (tT0SJout  - 8.8) 

(")  (PH)„U,  ■ 70 

(l)  [AlKlou<  • 2(CT>out  - |ACYlout 

(J)  [W)out  - lWlin  * [NaOH]addad  • 23  x 103 

(K)  No  Change:  Same  as  VI.K 

VIII.  Caustic  Neutral izat ion  (NaOH)/Using  H^O^,  ’ or  HN°3 

(A)  Calculate  (C^)  jn 

(1)  Analytical ly  determine  [ALK] [n  and  (PM)|n 

(2)  Calculate  (CT)  in  Using  Eqn.  II. A 
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(B)  Stoichiometric  Equation 

(1)  2NaOH  + H2SOi4  ->•  NaSO^  + 2H20 

(2)  NaOH  + HC1  -*•  NaCl  + H20 

(3)  NaOH  + HN03  -*•  NaN03  to  HjO 

(C)  Ca lculate  ly°t1added.  (lHCli,jded  or  [H"°,Iadded) 

(1)  During  neutralization  with  either  H2S0^,  HC1,  or 

HNV  <CT>out  ' <CT>in 

(2)  Assuming  neutralization  reaction  is  carried  to  pH  - 7 

( [ H ] = ,0  7),  then  using  Eqn.  II. A,  [H+]  » 10  , and 


(CT)  , calculate  [ALK] 
r out  out 

(3)  If  neutralizer  is  H2S0^,  then  [H2S0J)] 

[ALK]  )/2  (moles/ll) 
out 


([ALK] 


added 


(4)  If  neutralizer  is  either  HC1  or  HNO^  then  [HC ' 1 a<jded  ” 

( [ALkT.  - [ALK]  J,  [HNOj  . . . - f[ALK] . - [ALK]  ) 

in  out  3 added  in  out 

(D)  [SS]  „ = [SS],  (for  H.SO,  , HC1  6 HNOj 

out  in  L *4  3 

(E)  (pH)out  = 7,0  (for  H2S04 * HCI  6 HN°3 ) 

(F)  [ACY]  - 2(C_)  - (ALK)  (for  H SO, , HCI  6 HNO  ) 

out  T out  out  2 A 3 

(G)  [C0ND]Qut  - 1,786  ([TDS]out  - 8.8) 

(H)  If  H2S0^  is  neutralizing  agent 

">  'TDS'out  ■ ITDSlin*  'Wadd.d  * 98  * '°J 
<2>  ITSlou,  * [SS1out  * lT0Slout 

(3)  IS,Oout  * tS°i."'in  * ■W.dd.d  * 96  * lt|J 


BB MB 


(M 

No  Change: 

Sulfites  TNT 

Ni tr i tes 

TKN  TOC 

Chlorides  Ammonia 

Orthophosphates 

COO  MBAS 

Sod i urn 

Ni trates 

Ca 

if 

HC 1 is  neu tra 1 i zi ng 

agent 

(I) 

[TDS]  , ■= 
out 

lT0S1i„ 

♦ [HC1]  ^ * 

added 

36.5  x !03 

(2) 

[TS]  = 

out 

ISSlou. 

+ [TDS] 

out 

or  [TS] 

out 

- [TS],  + t HC 1 ] , 

in  added 

* 36.5  x 103 

(3) 

[Cl  ] « 

out 

[Cl  ] . 

in 

+ [HC 1 ] , * 

added 

35.5  x 103 

w 

No  Change: 

TKN 

Sod i urn 

Orthophosphate 

COD 

TNT 

Ca 

TOC 

Ammon i a 

MBAS 

Ni  trates 

Su 1 f i tes 

Ni tri tes 

Su I fates 

Color 

If 

HNO^  is  neut 

ra 1 i 2 i ng 

agent 

(1) 

[TDS]  - 

out 

[T0S]in 

♦ [HNVadd,d 

* 63  x 103 

(2) 

[TS]  „ » 
Jout 

tSS]out 

* |TDSlou, 

or  [TS!out 

- [TS]. 

„ * '""V.dd.d 

( * 63  x 103 

(3) 

[N°3>out  - 

[no3] 

* '""V.dd.d 

* 62  x 103 

(M 

No  Change: 

Sod i urn 

Nitrites 

TKN  MBAS 

TNT 

Color 

COD  Sulfite 

Ammon i a 

Orthophosphates 

TOC  Chloride 

Sul  fates 

Ca 

^3 


IX.  Caustic  Neutral i zat ion  (NaOH)/Using  CO^ 

(A)  Calculate  (CT). 

Tin 

Same  as  V I I I . A 

( B)  Stoichiometric  Equation 


2NaOH  + H2C0  -*■  Na2C0  + 2H20 

(C)  Calculate  [CO^^ 


(I)  During  neutralization  with  C02> 


(2) 


[ALK]  = [ALK], 
out  in 


Assuming  neutralization  reaction 
then  using  Eqn.  I I .A,  [H  ] = 10  ^ 


is  carried  to  pH  = 7, 

and  [ALK]  , calculate 
out 


(3)  ^Vadded  (<Vout  ' (Vin 


(D) 

rssU  = 

^in 

(E) 

(pH)out  = 

7-0 

(F) 

[ACY]  < 

out 

= 2 (CT)  „ - [ALK] 

T out  out 

(G) 

[COND] 

out 

= 1 .786  ( [TDS]  - 8.8) 

out 

(H) 

ITOS]  . 
out 

■ * gladded  . W x I0! 

(1) 

[TS]  = 

out 

tTDS1out  * ISS1ou, 

or  [TS]  , = [TS] . + [CO,]  ...  * 

out  in  2 added 


TKN 


TNT 


1 
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A. 2. 2.  SETTLING  MODEL 

The  SETTLE  model  takes  into  account  of  settling  ponds  or  clarifiers. 
A removal  factor  RF  is  computed  based  on  pond  area  and  flow  rate.  This 
removal  factor  is  applied  to  parameters  involving  settleable  solids 
according  to  the  following  formulas. 


| 

| 

I 
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Settle  Model 


(A)  Basic  Removal  Equation 

[CONC]  „ - [CONC] . 
out  in 


(1-0.82  exp  t-Q,/2?80A]  ) 
(1-0.00205  exp  [-Q,/2780A]  ) 


Removal  Factor 


or,  [CONC]  _ = [CONC].  * RF 
’ out  in 


where, 


Q.  = flow  in  million  gallons/day  (MGD) 

A = surface  area  of  clarifiers  (acres) 


(B)  Suspended  Sol  ids  (SS) 

[SS]  , = [SS].  * RF 

1 out  in 


(C)  Volatile  Suspended  Solids  (VSS) 

[VSS]  =>  [VSS],  * RF 
out  in 


(D)  Turbidi ty  (T) 

[T]  - [T] . * RF 

1 out  in 


(E)  Chemical  Oxygen  Demand  (C00) 

[COD]  „ - 0.3  [COD].  + 0.7  [COD].  * RF 
out  in  in 


(F)  Total  Organic  Carbon  (TOC) 

[T0Clout  * °'3  [T0Clin  + °’7  [T0C]in  * RF 


^7 


(G)  Total  Kjeldahl  Nitrogen  (TKN/N) 

[TKN/N ] = [NH  /N] . + ([TKN/N],  - [NH,/N] . ) * RF 

out  3 m in  3 in 

where , 

[NH  /N]  = ammonia  nitrogen 


(H)  Total  Solids  (TS) 

[TS]  = [TDS] . + [SS]  . * RF 

out  in  in 

where , 

[TDS]  = total  dissolved  solids 

(l)  Conduct i v i ty  (COND) 

[COND]  = 1.786  * ( [TDS]  . - 8.8) 

out  in 


(j ) No  Change: 

(1)  pH 

(2)  Total  Dissolved  Solids 

(3)  MBAS 

(A)  Sulfites 

(5)  Chlorides 

(6)  Sodium 

(7)  TNT 

(8)  Alka I in  i ty 

(9)  Acidity 

(10)  Ammonia 

(11)  N i trates 

(12)  Nitrates 

(13)  Color 

(14)  Orthophosphate 


A.  2.3  DOME  MODEL 


* 
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Domestic  or  sanitary  type  treatment  processes  are  simulated  by  the  DOME 
subroutine.  Two  major  types  of  domestic  treatment  systems  are  included  in  DOME. 
They  are  the  trickling  filter  and  activated  sludge  systems.  While  each  of  these 

systems  may  be  employed  in  many  different  mod i f i ca t ions  or  configurations, 

r 

DOME  is  equipped  to  handle  only  a few. 

DOME  is  suited  for  single  stage  filters  only.  Three  different  recycle 
configurations  can  be  simulated:  no  recycle;  recycle  from  the  filter  effluent; 

or  recycle  from  the  secondary  clarifier  effluent.  Recycle  conf igurat ion  is 
controlled  by  input  of  the  appropriate  values  of  R2,  R3,  and  R4;  the  trickling 
filter  recycle  terms  in  DOME.  Recycle  from  the  filter  effluent  is  R2;  recycle 
from  the  secondary  clarifier  is  R3;  and  recycle  to  the  filter  is  R4.  Obviously, 
if  RA  is  zero,  R2  and  R3  must  be  zero  and  DOME  will  simulate  a no  recycle  system. 
For  a system  with  recycle  from  the  filter  effluent,  R3  is  zero  and  R2  is  equal 
to  R k.  For  recycle  from  the  secondary  clarifier  effluent,  R2  is  zero  and  R3  is 
equal  to  R^. 

Activated  sludge  systems  simulated  by  DOME  include  the  conventional  or 
plug  flow  type  and  the  complete  mix  type  (which  includes  extended  aeration.) 

Con tact-s tab i 1 i za t ion , step-aeration,  or  other  such  modifications  of  the 
activated  sludge  process  are  not  simulated  by  DOME.  The  net  yield  term  supplied 
to  the  program  by  the  user  must  be  determined  for  the  actual  system  being 
simulated.  Calculation  of  the  net  yield  term  can  be  done  as  follows: 

Y = (Fw)  (5Sw)  + (F)  (SSo) 

n F (BOD . + BOD  ) 

i o 

Y = Net  yield  (lbs.  SS  produced/lb.  BOD  used) 
n 

F^  = Volume  of  sludge  wasted  (MGD) 

SS  = Suspended  solids  concentration  of  waste 
. W sludge  (mg/1) 


*♦9 


F * Plant  influent  flow  (MGD) 

SSq  * Plant  effluent  suspended  solids  (mg/1) 
BOD.  ■ Plant  Influent  BOD  (mg/1) 

B0Dq  ■ Plant  effluent  BOD  (mg/1) 


Information  necessary  for  this  calculation  can  be  obtained  from  the  plant's 
operating  logs  or  measured  during  the  preliminary  survey. 

Process  Equations  for  Subroutine  DOHE 
A listing  of  the  mathematical  expressions  used  in  the  formulation  of 
the  DOME  subroutine  follows: 

Primary  Clarifier 

RF  = 1.0  - 0. 82 EXP (-Q/2780A) 

1.0  - 0. 00205EXP (-Q/2780A) 


RF  = Removal  factor 
Q = Plant  flow  (MGD) 

A = Primary  clarifier  surface  area  (acres)  AREAP 


TSS 


out 


BOD 


out 


COD 


out 


TOC 


out 


VSS 

out 

Turb 

out 


TKN 


out 


TSS.  * RF 
1 n 

0.3BQD.  + 0.7B0D.  * RF 

in  in 

0.3C0D.  + 0.7C0D.  * RF 

in  in 

0.3T0C.  + 0.7T0C.  * RF 

in  in 

VSS.  * RF 
1 n 

1 Turb.  * RF 
1 n 

Ammonia,  + (TKN.  - Ammonia 
in  in 


!„>  * "F 


Assumes  BOD,  COD  and 
TOC  are  70%  particulate 
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Trickling  F i 1 ter 
BOD  reduction 


E = 1.0 

1.0  + 0.0085  [^]i 

E = Decimal  BOD  removal 

W = BOD  loading  (lbs/day)  to  filter  (not  including  recycle) 

V = Volume  of  filter  (acre-ft) 

F = Recirculation  factor  = 1 . 0 + R R = recirculation  ratio 

1.0  + 0.1R 

= recycle  flow 
plant  flow 


A temperature  correction  has  been  incorporated  in  the  BOD  reduction. 
Et=E20  (1.038)T*20 

E.j.  * Filter  BOD  removal  at  temperature  T 
E^  = Filter  BOD  removal  at  20  °C,  as  above 
T = Temperature  of  wastewater  (°C) 


| 

I 


BOD  = BOD.  * Et 
out  in  T 

SS  = SS.  + lbs,  solids  produced 

out  In  8.3*  (plant  flow  + recycle  flow) 

Tbs.  solids  produced  = 0.2  (lbs  BOD.  - lbs  BOD  ) 

in  out 


Nitrification 

E - EXP(-KD/Qn) 

E = Decimal  reduction  in  NH.-N  concentration 

K = Reaction  constant  related  to  the  specific  surface  of  the  media 
D - Depth  of  filter  (ft) 

Q = Hydraulic  loading  rate  (MGAD) 

n = Constant  related  to  the  specific  surface  and  configuration  of 
packing  of  the  media 
= 2.36  1"  rock  med i a 
= 3- 80  2iM  rock  media 

K.0  « K at  20°  C • 0.23  1"  rock  media 
0.13  2JM  rock  media 

K - K20(1.07)T  where  T is  wastewater  temperature  (°C) 
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Ammon i a 


out 


Ammonia.  * E 
i n 


N i trate 


out 


Nitrate,  + Ammonia.  - Ammonia 


'in 


i n 


OJt 


TKN 


» TKN,  - Ammonia,  + Ammonia 


out 


. ■ • - 1 I ruiaiis/i  i i a • _ 

out  in  In 

Trickling  filter  secondary  clarifier 

The  formulation  for  the  trickling  filter  secondary  clarifier 
as  that  for  the  primary  clarifier.  The  exception  being  that  only 
solids  removal  is  considered. 

Aeration  basin 


BOD  reduction 


is  the  same 
suspended 


BOD  = BOD.  - 
out  in 

MLSS  = Mixed  liquor  suspended  solids 
Vol  = Aeration  basin  volume  (MG) 

Y - Net  yield  - lb.  SS  produced/lb.  BOD  used  in  aeration 
Fn  » Plant  flow  + recycle  flow  (MGO) 

Theta  - Sludge  retention  time  or  sludge  age 


(MLSS)  (Vol) 
(Yn)  (Theta)  (f) 


1.0 

Y (PL  I ) 
n 


where  Pll  ■ (BOD.  ) (Plant  flow) 
i n 

(MLSS) (Vol ) 


Ni tri f i cat  ion 


A double  "halving"  or  Balzano  type  iteration 
concentrat ions  of  ammonia  and  the  nitrifying  bacte 


is  performed  on  the  effluent 
ria,  N i trosomonas , until  the 


governing  system  equations  are  satisfied: 


X 


i 


S. 

S' 


K 


o 


Y 

K 


r 


X 

X 


s 

i 


D 


o 

t 


* WA 

Y<lt5  * s„> 

K S X D 
root 

“5?  +”5 


Ammon  la  In 
Ammonia  out 

Synthesis  rate  constant 
Yield  coefficient 
Saturation  coefficient 
Ni trosomonas  in 
Ni trosomonas  out 
Detention  time  (Vol/Flow) 
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A detailed  explanation  and  justification  of  assumptions  may  be  found 
in  "A  Generalized  Computer  Model  for  Steady-State  Performance  of  the 
Activated  Sludge  Process",  Water  Pollution  Control  Research  Report,  by 
Robert  Smith  and  Richard  G.  Eilers,  FWPCA  (Washington,  1969). 

Activated  sludge  secondary  clarifier 

SSF  , (SS.HF)  (3.785) 

(24)  (A) 

2 

SSF  = Solids  surface  feed  (kg  SS/m  /hr 

F = Plant  flow  plus  recycle  flow  (MGD) 

SS.  = Mixed  liquor  suspended  solids  (mg/1) 

SS1  = 4.5  + 8.6  (SSF) 
out 

A . 3-0  Use  of  Plant  Site  Model: 

The  plant  site  can  be  modeled  and  various  operating  conditions  simulated 
using  a card  deck  containing  the  following  subroutines: 

MAIN:  sets  the  system  dimensions  and  calls  appropriate  subroutines; 

PROCES, START,  TOP. 

PROCES:  defined  in  previous  discussion.  This  program  reads  in  process 

mnemonics  for  processes  present  at  the  particular  site  being  studied. 

Appropriate  process  data  is  called  out  of  process  library.  Process  data  can 
be  altered  to  better  represent  a particular  site. 

Process  Library:  described  in  previous  discussion. 

START:  Reads  in;  the  number  of  processes  (sources),  the  number  of  •] 

branches,  the  number  of  parameters  and  the  number  of  outfalls  in  this  site 
study.  Also  reads  in  the  numbers  of  the  particular  pollutants  in  the  current 
study  (NPLIST).  Reads  in  branch  names.  Branches  are  named  for  convenient 
interpretation  of  output  data.  Calls  CHK1  . 
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CHKI : Reads  in  the  number  of  non  competing  parameters  (used  in  RM) ; reads 

in  the  names  of  all  of  the  pollutants  in  the  master  list,  fcr  each  the 
number  of  measurement  methods  available  and  the  sampling  frequency.  This 
program  also  prints  out  the  input  information  available  at  this  point.  In 
particular  this  include  the  process  data  and  the  topology  matrix. 

TOP:  This  program  was  described  in  detail  in  previous  paragraphs  but 

essentially  it  does  the  following:  reads  in  treatment  data;  takes  process 

data  and  site  topology  data  and  calculates  the  flow  and  pollutant  parameters 
for  each  branch  of  the  system.  The  effects  of  treatments  are  included  in 
the  calculation;  calls  the  treatment  subroutine  NEUTRA,  SETTLE,  and  DOME; 
call  CHK2  to  output  results. 

NEUTRA:  This  subroutine  is  called  by  TOP  if  acid  or  caustic  neutralization 

is  required  in  any  branch. 

SETTLE:  This  subroutine  is  called  by  TOP  if  a settling  pond  or  clarifier  is 

required  in  any  branch. 

DOME:  This  subroutine  is  called  by  TOP  if  domestic  or  sanitary  treatment  is 

requi red . 

CHK2:  This  program  essentially  writes  out  the  results  of  the  TOP  calculations, 

i.e.,  the  system  simulation. 

FLAG:  This  subroutine  is  not  used  in  system  simulation  but  in  system  planning 

with  RM.  However,  it  is  called  by  TOP. 
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B.  Resource  Matching 

After  the  survey  planner  has  modeled  the  plant  site  and  gone  through 
several  simulations  for  plant  operating  levels  probable  during  the  survey 
period,  he  needs  to  determine:  at  which  of  the  branches  should  he  take 

his  samples;  for  which  parameter  should  each  sample  be  analyzed;  for  the 
proposed  number  of  samples  and  parameters  what  equipment,  time,  and  per- 
sonnel are  required  or  conversely  for  the  available  time,  equipment,  and 
personnel  how  many  samples  and  parameters  can  be  considered.  This  is  the 
survey  planning  operation. 

This  survey  planning  is  accomplished  using  the  "site  modeling  and 
simulation"  program  discussed  in  the  previous  section  of  this  report  and 
the  resource  matching  (RM)  program  and  its  subroutines. 

Data  is  prepared  for  the  RM  program  by  TOP,  LEVEL,  and  FLAG  sub- 
routines. As  described  in  the  previous  section,  the  TOP  subroutine 
determines  the  flow  and  the  parameter  value  for  each  pollutant  at  each 
possible  measure  point  in  the  system.  The  LEVEL  subroutine  rearranges 
TOP  output  data  properly  for  use  in  RM. 

B.1.0  LEVEL  Subroutine: 

RM  considers  measure  points  in  subsets  or  "levels".  These  "levels" 
are  defined  such  that  if  a pollutant,  say  sulfates,  is  measured  at  all 
points  in  a "level"  then  all  the  sulfates  in  the  system  have  been  detected. 
The  following  paragraphs  explain  the  level  designation. 

The  three  dimensional  space  Pjjk  is  required  to  designate  the  pollutants 
at  all  points  throughout  the  site.  Note  the  following  system  topology  and 
pollutant  Pi  illustrated  in  Figure  IX. 


I 
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OUTFALL 

Level  Designations 
Figure  IX 

Pollutant  sources  are  at  points  1,  2,  3,  4,  5.  Whatever  the  source, 
if  Pi  is  measured  at  the  outfall,  one  sample  point  detects  all  of  Pi. 
However,  if  Pi  is  measured  at  level  2,  two  sample  points  are  required  to 
detect  all  of  Pi*  Thus,  pollutant  P is  designated  as 

pijk 

where  i is  for  the  ith  pollutant 
j is  for  the  jth  level 

k is  for  the  measure  point  at  the  j_th  level. 

One  needs  to  consider  the  pollutant  vectors  in  subsets  to  insure  that 
all  of  a particular  pollutant  has  been  detected.  For  example,  for  pollutant 
Pi  the  sets  are 

Pill 

or 

Pi  2 1 + Pi  22 

or 


P1 3 1 + p132  + P1 3 3 
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In  reading  pollutants  into  a computer  for  the  matching  program  the 
pollutants  will  have  a name  and  number  designation,  a measurement  level, 
a measurement  point  and  an  approximate  measure  number  (concentration,  pH, 
other).  For  example,  if  sulfate  is  the  number  5 pollutant,  sulfate  at 
the  fourth  measure  point  on  the  third  level  may  be  P (5,  3,  4)  - 50  mg/1. 

The  resource  matching  program  is  supplied  with  the  predicted  plant 
site  state  from  the  site-topology  model.  That  is,  a predicted  number  is 
given  for  Pijk  for  each  parameter  i,  for  each  level  j,  and  for  each  point 
k at  each  level. 

For  the  sewer  network  defined  by  the  A (1,J)  matrix,  the  LEVEL 
subroutine  generates  the  sets  of  points  that  form  levels. 

B.2.0  FLAG  Subroutine: 

The  points  where  samples  are  to  be  collected  and  the  parameters  for 
which  the  samples  are  to  be  analyzed  are  determined  by  the  RM  logic  or  by 
the  survey  planner  using  the  FLAG  subroutine. 

The  RM  logic  assigns  the  outfall  as  the  measure  point  for  all  para- 
meters unless  it  is  forced  to  seek  a higher  level  because  of  low 
concentrations  at  the  outfall  or  unless  it  is  forced  to  seek  a higher 
level  because  of  violation  of  a constraint.  When  a constraint  is 
violated,  the  RM  logic  seeks  new  methods  and/or  new  levels.  At  whatever 
level  measurements  are  to  be  taken  for  a given  parameter,  all  points  at 
that  level  must  be  measured. 

However,  as  a result  of  system  modeling  and  simulation  of  expected 
plant  operating  conditions  the  survey  planner  may  wish  to  take  additional 
measurements  to  obtain  further  useful  information.  By  the  use  of  the  FLAG 
subroutine  the  survey  can  cause  RM  to  consider  (flag)  additional  measure 


points  for  each  parameter. 
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The  RM  logic  assigns  resources  to  these  flagged  points  before  making 
assignments  through  its  own  logic.  A flag  designates  a particular 
pollutant  to  be  measured  at  a particular  measure  point.  Three  types 
of  flagging  procedures  are  available: 

1.  Automatic  flagging 

2.  Manual  flagging  of  sources 

3.  Manual  flagging  of  branches 

For  automatic  flagging  the  program  computes  the  total  mass  of  any 
pollutant  discharged  at  the  outfall  and  the  total  mass  for  that  pollutant 
from  each  process  (source) . The  program  will  then  automatically  flag  any 
source  which  contributes  more  than  X (decimal  fraction)  to  the  total  at 
the  outfall.  The  source  is  flagged  for  this  parameter  only.  The  survey 
officer  specifies  the  flagging  level  X for  each  parameter. 

The  FLAG  subroutine  also  permits  manual  flagging  in  two  parts.  The 
first  part  allows  the  user  to  select  specific  processes  (sources)  and 
parameters  he  wishes  to  flag.  The  second  part  allows  specific  branches 
and  parameters  to  be  flagged.  In  both  cases  the  process  (branch)  numbers 
and  corresponding  parameter  numbers  are  read  in  by  FLAG.  FLAG  also  outputs 
the  names  of  the  parameters  and  their  flagging  levels,  and  the  names  of 
the  processes,  branches  and  parameters  selected  by  the  user. 

B.3.0  Other  RM  Data  Required 

The  resources  available  for  conducting  the  survey  are  in  the  form  of 
the  Measurements  Resources  Model,  Sampler  Data  and  Flow  Measurement  Data. 

The  measurement  resources  model  consists  of  lists  of  measurement 
methods  applicable  for  each  pollutant.  There  are  up  to  three  methods  listed 
for  each  pollutant  in  the  master  list.  Different  methods  for  a given 
pollutant  may  Involve  an  entirely  different  technique  (such  as  Specific 
Ion  Electrode  or  132C  Phenate  for  Ammonia  measurement)  or  may  involve  the 
same  technique  performed  by  a different  class  of  analyst. 
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For  each  method  for  each  pollutant,  defining  characteristics  are  tabulated 
in  matrix  form  for  ready  access  by  a computer  program.  These  defining 
characteristics  are:  method  name,  minimum  acceptable  concentration, 

equipment  or  preparation  process  nstne,  equipment  or  preparation  time  per 
sample,  van  space  required,  analysts  time  required  for  each  of  four 
classes  of  analyst,  cost  per  sample,  sample  volume  per  sample,  at  which 
lab  analysis  is  done.  The  measurement  Resources  Model  is  illustrated 
in  Figure  X. 

The  RM  program  is  also  supplied  with:  the  number  of  measure 

methods  available  for  each  parameter;  the  sampling  frequency  for  each 
parameter;  the  names  of  the  pieces  of  equipment  in  all  resource  models 
and  time  constraint  on  each;  constraints  other  than  on  pieces  of  equipment. 

B.4.0  Resource  Matching  Model  Logic: 

The  function  of  the  resource  matching  model  is  to  take  the  predicted 
plant  site  state,  the  available  resources  and  the  constraints  and  to 
determine  a feasible  survey  plan.  The  plan  will  include  measure  points 
selected  by  the  RM  logic  and  points  flagged  by  the  survey  planner. 

The  predicted  plant  site  state  is  supplied  by  the  TOP  subroutine 
(plant  site  model)  and  consists  of  flow  and  pollutant  parameter  values 
for  each  branch  (possible  measure  point)  of  the  system.  The  plant  site 
state  is  determined  for  a specific  set  of  processes  and  operating  levels 
chosen  by  the  survey  planner  as  conditions  one  might  expect  to  occur  during 
the  survey  period. 

The  available  resource  data  is  supplied  by  the  resource  models. 

For  each  pollutant  parameter  to  be  measured  on  the  survey,  there  are  one 
to  three  measurement  methods  defined  in  the  resource  model  part  of  the 
program.  This  model  lists  the  equipment,  personnel,  and  time  required 

for  a given  measurement.  I 

, J 
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Other  resource  data  was  discussed  in  the  previous  section  of  this  report 
and  is  illustrated  in  Figure  X. 

The  constraints  are  generally  in  the  form  of  limited  time  available 
for  use  of  a given  piece  of  equipment  or  class  of  analyst.  Limits  on 
available  space  and  cost  can  also  be  included. 

The  survey  plan  determined  by  the  resource  matching  model  includes: 
specification  of  a measurement  method  for  each  parameter;  specification  of 
the  measure  points  for  each  parameter;  specification  of  resources  used  by 
measurement  of  each  parameter;  a listing  of  the  total  commitment  for  the 
survey  period  of  each  piece  of  equipment  and  each  class  of  analyst;  a 
listing  of  the  constraint  violations. 

A feasible  survey  plan  is  feasible  in  the  sense  that  limits  on 
equipment  and  personnel  time  available  have  not  been  exceeded.  If  no 
feasible  alternatives  are  available,  the  program  will  specify  a non- 
feasible  survey  plan  and  note  the  constraints  violated. 

The  survey  planner  can  generate  alternative  feasible  and  non-feasible 
survey  plans  by  altering  priorities  in  measurement  methods,  priorities  in 
the  order  of  assigning  resources  to  parameters,  or  by  varying  constraints 
on  time,  equipment  or  personnel.  Also,  different  sets  of  measure  points 
may  be  flagged  by  the  survey  planner. 

Alternative  solutions  are  possible  because  different  measurement 
methods  have  different  requirements.  The  characteristics  of  the  measurement 

methods  which  afford  a basis  for  trade  off  are: 

1.  Different  measurement  methods  may  require  different  minimum 
pollutant  concentration. 

2.  Different  measurement  methods  may  require  different  amounts  of 
time  expenditure  by  technicians  and  by  equipment. 


3.  Different  methods  require  different  pieces  of  equipment,  and 

there  are  time  and  space  constraints  on  personnel  and  equipment. 
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The  programmer's  manual  in  the  appendix  includes  a flow  chart  or  block 
diagram  of  the  re>'  mrce  matching  logic.  The  following  paragraphs  describe 
that  logic  in  general  terms. 

The  resource  matching  program  is  supplied  with  a predicted  plant 
site  state  from  the  site-topology  model  (TOP).  That  is,  a predicted 
number  is  given  for  for  each  parameter  i,  for  each  level  j,  for 

each  point  k at  that  level.  Also,  for  each  parameter  to  be  measured  on 

I 

a survey,  on  to  three  measure  methods  are  defined  by  a measure  matrix 
pimi*  That  is,  methods  Ml,  M2,  M3  for  each  pollutant  i.  As  part  of 
the  method  definition,  the  minimum  parameter  value  is  given.  For  any 
lower  number,  the  method  loses  credibility. 

The  first  step  in  the  resource  matching  logic  is  the  generation  of  a 
set  of  matrices  which  will  list  all  feasible  measurement  method  matrices 
(feasible  in  terms  of  minimum  parameter  value)  for  each  P^^.  Every 
number  P^jk  is  compared  with  the  minimum  parameter  value  in  each 
appropriate  measure  method  matrix  P^m£.  If  the  pijk  nurober  is  equal  to 
or  greater  than  the  required  minimum,  the  P^jj  is  saved  or  stored  as 
Pijkj.  This  is  a feasible  measure  method  matrix  for  parameter  i,  level 
j,  point  k,  method  i.  Each  matrix  is  of  the  form  of  the  previously 
defined  method  matrix  P^mj  illustrated  in  Figure  X.  The  complete  set 
of  matrices  contains  all  feasible  methods  for  each  parameter  at  each 
measure  point  in  the  system.  The  resource  matching  model  will  next 
determine  a feasible  survey  schedule.  The  schedule  is  the  assignment  of 
measure  point  or  points  for  each  parameter  and  a measure  method  for 
each  parameter  at  those  points.  The  schedule  is  feasible  if  the  total 
commitment  of  time  and  equipment  does  not  exceed  availability  limits. 

It  is  reasonable  to  assume  that  for  any  system  there  are  many  feasible 
surveillance  schedules. 
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The  survey  officer  may  choose  a best  schedule  in  terms  of:  value  in  dollars, 

technician  time,  space  required,  time  used  on  automated  equipment,  or  total 
survey  time. 

The  feasible  solutions  are  determined  by  assigning  a measurement  method 
to  each  pollutant  at  the  outfall  or  at  the  lowest  level  at  which  a Pijk£ 
exists  for  that  pollutant.  This  is  started  with  some  arbitrary  PI.  As 
each  additional  pollutant  is  added  an  account  of  the  commitment  of 
resources  is  kept.  The  limits  or  constraints  on  the  resources  are  known. 

When  any  constraint  is  encountered  the  program  first  seeks  an  alternate 
measurement  method  at  the  same  topological  level.  If  the  alternate  method 
is  found  and  does  not  violate  the  same  or  any  new  constraints  the  next 
pollutant  is  considered.  If  no  alternative  measurement  method  is  found 
at  the  same  level,  then  a new  level  is  tried.  It  no  alternative  is  found 
at  new  levels,  then  the  previously  considered  pollutant  is  tried  for  an 
alternative  method.  If  the  alternative  exists  Pj_l  is  replaced  by  P^  and 
the  next  pollutant  is  considered.  If  Pj_i  does  not  afford  a solution 
then  P-j_2  considered  and  so  forth.  If  no  feasible  alternative  exists 
the  program  returns  to  the  Pj  entered  with  the  violation.  The  assignment 
continues  with  the  violation  noted.  Alternative  solutions  are  generated 
by  selecting  alternative  sequences  for  considering  measurement  methods 
or  for  considering  the  different  pollutant  parameters. 

The  systems  analysis  of  the  survey  design  will  reduce  the  time 
required  for  survey  planning.  Also,  the  systems  approach  will  allow 
consideration  of  many  alternative  plans,  with  each  of  these  plans  more 
precisely  defined  quantitatively.  It  is  reasonable  to  expect  that  a 
computerized  system  analysis  will  be  used  to  seek  alternatives  during  the 
survey,  as  well  as  during  planning. 
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B.5.0  Output  Data: 

Some  of  the  output  of  the  RM  program  is  a reiteration  of  the  input 
data.  Other  output  data  is  essentially  the  survey  design  and  is  discussed 
in  the  following  paragraphs.  Sample  output  data  is  given  in  the  appendix 
under  Program  Listing  and  in  the  User’s  manual. 

Each  branch  of  the  system  is  listed  and  along  with  each  parameter 
that  is  to  be  measured  at  that  branch.  Also  for  each  branch  the  sampling 
procedure  is  defined  in  terms  of  the  number  of  samples  taken  at  that 
point,  the  sampling  method,  and  the  total  technician  time  required  for 
those  samples. 

Next,  for  each  parameter  in  the  study,  the  points  where  this  parameter 
is  to  be  measured  are  listed.  At  each  point,  the  measurement  method  and 
the  expected  value  are  given  along  with  the  total  number  of  samples. 

Items  of  equipment  used,  total  time  used  for  this  parameter,  and  any 
equipment  time  constraint  are  noted  for  each  parameter.  In  addition,  the 
total  time  commitment  for  this  parameter  for  each  of  four  classes  of 
analysts  is  listed  along  with  time  constraints  encountered. 

Next  the  survey  commitments  are  summarized  in  terms  of  the  demand  on 
each  piece  of  equipment  and  each  class  of  analyst. 

For  each  item  of  equipment  used  the  program  lists:  the  number  of 

parameters  using  this  item;  the  total  time  used;  was  constraint  violated; 
the  amount  of  the  violation;  space  required. 

For  each  class  of  analyst  the  program  lists:  the  total  time  commitment; 

was  constraint  violated;  the  amount  of  the  violation. 

Data  on  van  space  utilization  and  cost  can  be  read  out  if  proper 
data  is  input. 


B.6.0.  Subroutines  Required  for  RM 


MAIN  and  all  of  the  subroutines  required  in  Site  Modeling  and 
Simulation  are  required  in  Resource  Matching  (RM) . In  addition  the 
following  subroutines  are  required. 

LEVEL  arranges  TOP  output  for  use  in  RM. 

FLAG  allows  survey  planner  to  specify  for  RM  measure  points  in 
addition  to  those  selected  by  RM  logic. 

RM  reads  in  and  prints  out  resource  information,  allocates  resources 
and  calls  subroutines  when  needed  to  relieve  violations,  and  prints 
out  resource  allocation  information. 

CONCK  determines  feasible  measure  methods  and  allocates  resources  to 
flagged  points. 

TO  zeros  out  temporary  allocations  for  equipment  times. 

EQCHEC  determines  if  any  equipment  time  violations  occur. 

VSCHEC  checks  for  van  space  constraint  violations. 

CNCHEC  checks  for  analyst  time  constraint  violation. 

ADD  adds  current  allocations  of  equipment  time,  analysts'  times, 
van  space  and  costs. 

SET  stores  the  array  of  methods  used  at  a branch  in  a level  for  a 
given  parameter. 

PIEQCH  searches  past  assignments  to  attempt  to  relieve  constraint 
violations. 

SUBT  subtracts  current  allocations. 

INFORM  prints  out  branch  names  and  parameters  to  be  measured  at 
that  branch. 

SAMPLE  allocates  analyst  type  1 time  required  for  taking  samples  and 
measuring  flow. 

PRPAR  outputs  parameter  measurement  information. 
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USECT  totals  usage  equipment. 

TICHAN  converts  minutes  to  hours  and  minutes. 


66 


C.  MODEL  VERIFICATION  AND/OR  UPDATING 
(Subroutine  CORRCT) 

C.1.0.  INTRODUCTION 

The  CORRCT  Subroutine  accepts  site  (Measurement  data,  compares 
it  with  model  predicted  measurement  data  and  computes  and  prints  out 
suggested  new  values  for  the  parameters  of  the  process  models. 

The  subroutine  CORRCT  is  used  in  conjunction  with  the  program 
subroutines  which  make  up  the  site  model.  They  are  MAIN,  PROCES, 

START,  TOP,  NEUTRA,  SETTLE,  and  DOME. 

The  site  model  package  with  CORRCT  has  two  objectives.  First,  the 
process  models  can  be  updated  or  imporved  using  preliminary  survey  data. 

While  the  best  data  available  were  used  in  constructing  the  process  mo- 

1 

dels,  there  was  considerable  variation  in  data  for  a given  process  from 
one  plant  site  to  another  or  from  one  season  to  another.  Better  process 
model  values  may  be  obtained  with  CORRCT  using  measurements  from  judi- 
ciously selected  grab  samples  taken  on  a preliminary  survey.  The  out- 
put of  CORRCT  lists  new  values  for  certain  process  parameters.  If  the 
survey  officer  wishes,  he  may  enter  these  new  values  into  the  process 
models  using  procedures  outlined  in  the  discussion  of  the  PROCES  sub- 
routine . 

A second  use  of  site  model  package  with  CORRCT  would  be  verifica- 
tion of  the  model.  After  the  survey,  measurements  are  available  for 
obtaining  max,  min,  average  or  other  expected  values  for  each  para- 
meter at  each  branch  for  known  levels  of  production.  These  values 
can  be  entered  into  the  site  model  package  with  CORRCT  to  determine 
values  for  process  model  parameters  which  can  be  used  to  verify  the 


r 


67 


model  or  refine  the  model  as  appropriate.  With  confidence  in  the  site 
model,  one  can  now  use  it  to  predict  levels  of  potential  pollutants  for 
production  levels  other  than  those  which  occurred  during  the  survey. 

This  should  give  added  significance  to  the  survey  report. 

Ideally,  if  measurements  are  to  be  made  to  verify  the  process  mo- 
dels, the  measure  points  should  be  chosen  such  that  each  process  will 
be  determined  uniquely  by  the  measurements.  However,  if  the  survey 
officer  who  determines  the  measure  point  does  not  visualize  the  system 
in  terms  of  the  same  processes  as  are  used  in  the  process  model,  mea- 
surements may  not  uniquely  determine  process  parameters.  The  program 
is  designed  to  accept  such  non-unique  measurement  data  and  uniquely 
determine  processes,  where  possible,  and  make  approximations  on  other 
processes  informing  the  user  if  the  determination  is  non-unique. 

The  measure  points  for  uniquely  determining  all  processes  can  be 
identified  by  referring  to  the  linear  graph  of  the  system  (Figure  1 
for  example)  which  displays  all  processes  and  all  possible  measure 
points.  On  the  linear  graph,  identify  all  of  the  vertices  (junctions) 
mearest  to  the  processes.  If  at  a vertex  there  are  N branches  in- 
cident to  a vertex,  then  measurements  must  be  taken  or  implied  at 
N - 1 of  those  branches.  By  taking  N - 1 measurements,  the  Nth  mea- 
surement is  implied.  Note  that  one  of  the  N - 1 measurements  at  ver- 
tex A may  be  implied  by  measurements  at  vertex  B or  other  vertices. 

C . 2 . 0 . INPUT  DATA  AND  OUTPUT  DATA 

After  the  decisions  are  made  which  determine  where  measurements  j 

are  to  be  made,  this  information  is  supplied  to  the  CORRCT  subroutine 


in  the  form  of  the  MESUR(I)  matrix.  The  MEFr  (I)  matrix  contains  a 
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"0"  for  each  branch  where  measurements  are  not  taken  and  a "1"  for 
each  branch  where  measurements  are  taken. 

It  may  not  be  appropriate  or  even  possible  to  check  every  para- 
meter in  a process  model.  The  subroutine  CORRCT  only  applies  to  con- 
servative parameter  (i.e.  only  those  measured  in  terms  of  mass,  mg/i). 
The  program  allows  the  user  to  select  the  one,  two  or  several  (up  to 
20)  parameters  to  be  included  in  each  model  verification.  The  para- 
meters to  be  considered  in  a particular  study  are  Identified  by  their 
numbers  in  the  master  list  of  Table  II.  The  flows  must  be  included 
in  any  branch  where  other  parameters  are  considered. 

The  number  of  parameters  and  the  parameter  numbers  are  supplied 
to  CORRCT  for  each  branch  where  measurements  are  taken  using  BPLIST 
(IB,N),  where  IB  is  the  branch  number  and  N is  the  parameter  number. 
The  first  number  N for  each  branch  is  the  number  of  parameters  con- 
sidered in  that  branch.  The  subsequent  numbers  N are  the  parameter 
numbers  for  that  branch. 

The  next  inputs  to  CORRCT  are  the  measured  parameter  values 
YM(IB,N).  For  each  branch  IB  the  first  value,  for  N = 1,  must  be 
the  flow.  The  remaining  YM(IB,N)  ir-put  for  a branch  are  the  measured 
values  of  the  parameters  listed  in  .IPLIST(IB.N)  . 

The  subroutine  CORRCT  output  1 lsts  each  source  (process)  by 
number,  notes  if  the  determination  is  unique  or  estimated , and  gives 
the  old  parameter  and  the  new  parameter  for  each  parameter  verified. 

If  the  determination  is  unique,  the  branch  where  the  defining  mea- 


surement was  taken  is  listed. 
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C. 3.0.  SUBROUTINE  LOGIC 

The  logic  of  the  program  is  outlined  in  this  and  the  following 
paragraphs.  The  program  scans  each  branch  of  the  system  starting  with 
branch  number  one.  It  determines  how  many  sources  (processes)  contri- 
bute to  this  branch  and  specifically  which  sources.  The  program  also 
notes  if  any  of  the  contributing  sources  were  treated  on  or  before 
reaching  their  branch.  If  only  one  source  contributes  to  this  branch 
and  it  is  not  treated,  then  the  branch  measurements  are  taken  as  the 
suggested  new  values  for  the  source  (process)  parameters.  If  only 
one  source  contributes  to  this  branch  and  it  is  treated,  then  the 
effect  of  the  treatment  on  each  parameter  is  determined  using  the  ori- 
ginal model.  The  branch  measurements  are  modified  by  these  factors  to 
give  the  suggested  new  values  for  the  source  (process)  parameters. 

If  more  than  one  source  contributes  to  a branch,  the  program  de- 
termines if  all  but  one  has  previously  determined  uniquely.  If  all 
but  one  has,  then  the  suggested  new  parameter  values  are  obtained  by 
subtracting  known  source  values  from  the  measurements  to  give  a new 
uniquely  determined  source.  These  subtractions  are  done  in  terms  of 
mass  and  not  concentrations. 

If  more  than  one  source  contributes  to  a branch  and  more  than  one 
of  these  sources  has  not  been  previously  uniquely  determined,  then 
these  undetermined  sources  can  only  be  estimated.  For  each  undetermined 
source  J contributing  to  branch  I the  percent  contribution  for  parameter 
K is  determined  from  the  original  model.  The  branch  measurements  are 
modified  by  this  % to  give  estimated  new  values  for  the  source  parameter 
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K.  These  percentages  are  determined  in  terms  of  mass  and  not  concentrations. 

Due  to  the  treatment  there  may  be  some  nonlineraties  and  therefore 
some  errors  in  the  uniquely  determined  sources.  If  refinement  is  desired, 
the  suggested  new  values  could  be  assigned  to  the  processes  and  in  several 
iterations  these  errors  could  be  eliminated. 

Errors  in  the  estimates  for  the  non  unique  sources  could  have  many 
causes.  Attempted  improvements  through  interative  solutions  may  not  be 
fruitful.  If  better  results  are  desired,  the  best  course  would  be  to  re- 
assess the  selection  of  measurement  points. 
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I I I INDICATOR  MODEL 
A.  INTRODUCTION 

The  indicator  model  is  used  to  calculate  certain  parameters,  generally 
termed  process  variables  or  operating  character i s t i cs , that  indicate  the 
operational  state  of  the  system.  These  calculated  values  can  then  be 
compared  with  the  desirable  values  as  cited  in  design  books  or  manuals*. 

A comparison  will  give  the  engineer  insight  into  the  operation  of  the 
system  such  that  problem  areas  and  possible  causes  of  these  problems  are 
more  readily  apparent.  It  can  also  indicate  whether  inadequate  plant 
performance  is  due  to  design  deficiencies  or  poor  operational  practices. 

An  example  should  suffice  to  illustrate  this  more  clearly.  If  a 
plant's  effluent  suspended  solids  concentration  is  too  high,  the  first 
thing  the  engineer  might  suspect  is  the  secondary  clarifier.  Solids 
removal  efficiency  calculated  by  the  model  will  indicate  if  the  secondary 
clarifier  is  performing  poorly.  If  the  secondary  does  appear  to  be 
performing  poorly,  operating  characteristics  such  as  surface  settling 
rate  and  solids  loading  rate,  calculated  by  the  model,  can  be  compared 
to  those  suggested  in  the  leterature  as  above.  If  the  clarifier  is 
failing  due  to  a hydraulic  overload,  the  comparison  of  the  surface  set- 
tling rare  to  suggested  values  should  indicate  this.  If  the  operational 
charac ter  i st ics  of  the  secondary  clarifier  compare  favorably  with  those 
in  the  literature,  but  it  still  appears  to  be  failing,  an  investigation 
of  the  operational  characteristics  of  the  other  unit  processes  might  be 

Wastewater  Engineering,  Metcalf  and  Eddy,  Inc.,  McGraw  Hill  Book  Company, 

N.  Y.,  M.  Y.,  (1972) 

Process  Design  Manual  for  Upgrading  Existing  Wastewater  Treatment  Plant, 

U.  S.  EPA,  Office  of  Technology  Transfer,  (Washington  IS7M 


Environmental  Systems  Engineering,  L.  G.  Ric^,  McGraw  Hill  Book  Company, 

N.  Y.  N.'  Y.  , (1973) 
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helpful.  Should  the  process  loading  intensity  fo  the  aeration  basin  not 
compare  favorably  with  cited  values,  a poorly  settling  sludge  might  be 
the  cause  of  the  secondary  clarifier's  failure. 
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A. 1.0  OPERATIONAL  CHARACTERISTICS 

Like  the  treatment  system  model,  the  indicator  model  is  applicable 
to  trickling  filter  and  activated  sludge  type  treatment  plants.  Opera- 
tional  characteristics  calculated  by  the  model  are  as  follows: 

Primary  Clarifier 
Detention  time  (hrs.) 

Surface  settling  rate  (GPD/ft2) 

Weir  overflow  rate  (GPD/ft) 

Solids  loading  rate  (lb  SS/ft2  -day) 

Solids  removal  efficiency  (%) 

Trickling  Filter 
Hydraulic  loading  (MGAD) 

Organic  loading  (lb  BOD/day/acre-f t) 

Recirculation  ratio 

BOD  removal  efficiency  (%) 

Aeration  Basin 
Detention  time  (hrs) 

Volumetric  loading  (lb  BOD/day/lOOO  ft3) 

Process  loading  intensity  (lb  BOD/day/lb  MLSS) 

Sludge  age  or  retention  time  (days) 

Diffused  aeration 

Volumetric  air  supply  (CFM/1000  ft3) 

Organic  air  supply  (ft3/lb  BOD) 

Mechanical  aeration 

Volumetric  power  supply  (HP/1000  ft3) 

Organic  oxygen  supply  (lb  oxygen/lb  BOD) 


BOD  removal  efficienty  (%) 
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Secondary  Clarifier 
Detention  time  (hrs) 

Surface  settling  rate  (GPD/ft2) 

Weir  overflow  rate  (GPD/ft) 

Solids  loading  rate  (lb  SS/day/ft2) 

Solids  removal  efficienty  (%) 

Overal  Plant 

BOD  removal  efficienty  (%) 

Solids  removal  efficienty  (%) 

A. 2.0  MATHEMATICAL  RELATIONSHIPS 

Each  of  the  operational  characteristics  calculated  by  the  model  is 
represented  by  a mathematical  expression.  Many  of  these  expressions  are 
evaluated  at  both  average  and  peak  flow  conditions.  This  is  done  by  using 
either  the  average  of  peak  flow  term  in  the  expression.  However,  this 
may  not  always  be  completely  accurate,  as  the  BOD  and  suspended  solids 
concentrations  also  vary.  Such  error  may  be  overcome,  if  desired,  by 

substituting  the  peak  values  of  these  parameters  (BOD,  flow,  etc.)  for 

the  average  values  and  running  the  model  a second  time.  Average  values 
of  the  operational  characteristics  calculated  on  the  second  run  will  then 
represent  the  peak  values.  Following  is  a listing  of  the  relationships 

Primary  Clarifier 

Detention  time:  DT  *=  V 

1000  F/24 

V * clarifier  volume  (1000  gal) 

F = plant  flow  (MGD) 

Surface  settling  rate: 

SSR  = F * 106 
A 

F * plant  flow  (MGD) 

A = clarifier  surface  area  (ft2) 
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Weir  overflow  rate: 

WOR  = F * 106 
L 

F - plant  flow  (MGD) 
L = weir  length  (ft) 

Solids  loading: 

SL  ■=  (SS)  (F)  (8.34) 
A 


SS  » influent  suspended  solids  (mg/1) 
F - plant  flow  (MGD) 

A = clarifier  surface  area  (ft2) 


Solids  removal  efficiency: 


SRE 


SSI  - SSO 
SSI 


* 100. 


SSI  = plant  influent  suspended  solids  (mg/1) 
SSO  = primary  effluent  suspended  solids  (mg/1) 


Trickling  Filter 

Hydraulic  loading: 

HL  - F + R 
A 

F = plant  flow  (MGD) 

R = recycle  flow  to  filter  (MGD) 
A **  filter  surface  area  (acres) 

Organic  loading: 

OL  - (F)  (BOD)  (8.34) 

V 

F « plant  flow  (MGD) 

BOD  « primary  effluent  BOD  (mg/1) 
V » filter  volume  (acre-ft) 

Recirculation  ratio: 

RR  - R 
F 


F “ plant  flow  (MGD) 

R ■ recycle  flow  to  filter  (MGD) 
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BOD  removal  efficiency: 

E - BOD*  - BOD0  * 100 
BODj^ 

BOD^  **  primary  effluent  BOD  (mg/1) 
BODQ  » plant  effluent  BOD  (mg/1) 

Aeration  Basin 


Detention  time: 

DT  “ 24V 

(F+R) 1000 

V - aeration  basin  volume  (1000  gal) 

F * plant  flow  (MGD) 

R * recycle  flow  (MGD) 

Volumetric  loading: 

VL  - (BOD)  (F)  (8.34)  (7.481) 

V 

BOD  » primary  effluent  BOD  (mg/1) 

F « plant  flow  (MGD) 

V - aeration  basin  volume  (1000  gal) 

Process  loading  intensity: 

PLI  - (BOD)  (F)  (8.34)  (1000) 

(MLSS)  (V)  (8.34) 

BOD  - primary  effluent  BID  (mg/1) 

F • plant  flow  (MGD) 

MLSS  “ mixed  liquor  suspended  solids  (mg/1) 

V • aeration  basin  volume  (1000  gal) 

Sludge  age: 

SA  - (MISS)  (0.001V)  (8.34) 

(Fw)  (RSS)  (8.34)  + (F-Fw)  (SS)  (8.34) 

MLSS  » mixed  liquor  suspended  solids  (mg/1) 

V - aeration  basin  volume  (1000  gal) 

Fw  ■ waste  sludge  flow  rate  (MGD) 

RSS  ■ return  sludge  suspended  solids  (mg/1) 

F - plant  flow  (MGD) 

SS  ■ plant  effluent  suspended  solids  (mg/1) 

A basic  assumption  made  in  the  sludge  age  formulation  is  that  sludge 
wasting  is  from  the  sludge  recycle  line.  If  the  wasting  is  from  the 
aeration  basin,  the  actual  waste  sludge  flow  rate  must  be  “djusted 
before  input  to  the  model. 


(F  ) (MLSS) 

F = — - 

i RSS 

F^  = input  waste  sludge  flow  (GPD) 

Fa  = actual  waste  sludge  flow  (GPD) 

MLSS  = mixed  liquor  suspended  solids  (mg/1) 
RSS  = return  sludge  suspended  solids  (mg/1) 


Volumetric  air  supply: 


_ (AV)  (1000)  (7.481) 

A ‘ (24)  (60)  (V) 

3 

AV  = volume  of  air  delivered  (1000  ft  /day) 
V = aeration  basin  volume  (1000  gal) 


Organic  air  supply: 


OAS 


(AV)  (1000) 
(BOD)  (F)  (8.34) 


3 , 

AV  = volume  of  air  delivered  (1000  ft  /day) 
BOD  = primary  effluent  BOD  (mg/1) 

F = plant  flow  (MGD) 


Volumetric  power  supply: 


VPS 


7.481  HP 
V 


HP  = total  aerator  horsepower  (HP) 

V = aeration  basin  volume  (1000  gal) 


Organic  oxygen  supply: 

_ (N)  (HP)  (HR) 

(BOD)  (F)  (8.34) 

N = oxygen  transfer  rate  (lb  02/HP-hr) 
HP  = total  aerator  horsepower  (HP) 

HR  = hours  of  aerator  operation  (hr/day) 
BOD  = primary  effluent  BOD  (mg/1) 

F = plant  flow  (MGD) 

N = (2.5)  (1.024)  1-20  (.08) 


CS  = D.0.  saturation  in  aeration  basin  (mg/1) 
C = D.0.  cone,  in  aeration  basin  (mg/L) 

T = plant  effluent  temperature  (°C) 


A basic  assumption  made  in  the  organic  oxygen  supply  formulation  is 
that  the  mechanical  aerators  employed  are  rated  at  2.5  lb  02/HP-hr 
under  standard  conditions  as  defined  in  [1].  Also  assumed  is  an 
alpha  value  for  the  wastewater  of  0.8.  If  the  engineer  has  the 
actual  values  of  these  parameters  for  the  system  being  investigated, 
the  model  should  be  changed  accordingly. 


[1]  Wastewater  Engineering,  Metcalf  and  Eddy,  Inc.,  McGraw-Hill 
Book  Company,  N.Y.  , N.Y.,  (1972). 
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BOD  removal  efficiency: 

E * BOD*  - B0Do  *100. 

BOD^ 

BOD^  * primary  effluent  BOD  (mg/1) 

BOD0  - plant  effluent  BOD  (mg/1) 

Secondary  Clarifier 

For  the  secondary  clarifier,  the  flow  to  the  clarifier  depends  on 
whether  it  is  part  of  a trickling  filter  or  activated  sludge  type 
system.  For  the  trickling  filter  system,  the  flow  to  the  clarifier 
is  equal  to  the  plant  flow  plus  the  recycle  to  the  filter  minus  the 
recycle  from  the  filter  effluent.  For  the  activated  sludge  system, 
the  flow  to  the  clarifier  is  equal  to  the  plant  flow  plus  the  recycle 
sludge  flow. 

Detention  time: 

DT  - V 

1000  F/24 

V - clarifier  volume  (1000  gal) 

F - flow  to  clarifier  (MGD) 

Surface  settling  rate: 

SSR  - F * 106 
A 

F “ flow  to  clarifier  (MGD) 

A - clarifier  surface  area  (ft2) 

Weir  overflow  rate: 

WOR  - F * 106 
L 

F ■ flow  to  clarifier  (MGD) 

L - weir  length  (ft) 

Solids  loading  rate: 

SL  - (SS)  (F)  (8.34) 

A 

SS  - aeration  basin  or  trickling  filter 
effluent  suspended  solids  (mg/1) 

F • flow  to  clarifier  (MGD) 


d. 


Solids  removal  efficiency: 


79 


SSI  - sso 
sso 


*100 


SSI  - aeration  basin  or  trickling  filter 
effluent  suspended  solids  (mg/1) 

SSO  - plant  effluent  suspended  solids  (mg/1) 

A. 3.0  Data  Input 

The  indicator  model  is  designed  to  use  data  that  is  available  prior 
to  a survey's  implementation  and  also  the  data  generated  by  the  survey. 
Design  data  needed  should,  of  course,  always  be  available.  Operational 
data  required  will  be  available  from  the  daily  and  monthly  operation 
logs  as  required  by  Technical  Manual  5-665,  and  then  from  the  data 
generated  by  the  survey.  (If  the  value  of  any  parameter  is  not  available, 
it  should  be  replaced  by  -1.0  in  the  appropriate  columns  for  the  data 
input.)  Following  is  a listing  of  the  data  input  necessary  to  implement 
the  indicator  model. 


Trickling  Filter  Type  Plant 
Card  1 


Column  1:  Plant  type;  1 for  trickling  filter  type  plant. 

Columns  2-41:  Name  of  installation  e.g.  Fort  Ord,  CA 

(Should  be  centered) 

Card  2 - Primary  Clarifier 

Columns  1-8:  Total  volume  (1000  gal) 

(2  clarifiers  @ 4000  gal  ea.  - 8.0) 

Columns  9-16:  Total  surface  area  (Ft2) 

(2  clarifiers  0 200  Ft2  ea.  - 400.) 

Columns  17-24:  Total  weir  length  (Ft) 

(2  clarifiers  @ 500  Ft  ea.  - 1000.) 

Columns  25-32:  Design  flow  rate  (MGD) 

Columns  33-40:  Design  detention  time  (Hrs) 

If  not  readily  available,  calculate 

DT  - Total  Vol  (1000  gal)  X 0.024 


Design  flow  (MGD) 
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Columns  41-48: 


Columns  49-56: 


Design  surface  settling  rate  (GPD/Ft2) 
If  not  readily  available,  calculate 

SSR  « Design  flow  (MGD)  X 106 
Total  surface  area  (Ft2) 

Design  weir  overflow  rate  (GPD/Ft) 

If  not  readily  available,  calculate 

WOR  = Design  Flow  (MGD)  X 106 
Total  weir  length  (Ft) 


Card  3 - Trickling  Filter 


Columns  1-8:  Design  flow  rate  (MGD) 


Columns  9-16:  Design  organic  loading  (Lb  BOD/Acre-Ft  Day) 


Columns  17-24: 


Total  volume  (Acre-Ft) 

(2  filters  @ 3 acre-ft  ea.  • 6.) 


Columns  25-32: 


Total  area  (Acres) 

(2  filters  @ .75  acres  ea.  *»  1.5) 


Columns  33-40: 


Total  design  recycle  per  day  (MG) 

If  not  readily  available,  calculate 

R = Design  percent  recycle  X 100.  X Design  flow  (MGD) 


Card  4 - Secondary  Clarifier 


Columns  1-8: 


Columns  9-16: 


Total  volume  (1000  Gal) 

(2  clarifiers  @ 4000  gal  ea.  * 8.0) 

Total  surface  area  (Ft2) 

(2  clarifiers  @ 200  Ft2  ea.  » 400.) 


Columns  17-24: 


Total  weir  length  (Ft) 

(2  clarifiers  @ 500  Ft  ea.  - 1000.) 


Columns  25-32:  Design  flow  rate  (MGD) 


Columns  33-40: 


Columns  41-48: 


Design  detention  time  (Hrs) 

If  not  readily  available,  calculate 

DT  - Total  Vol  (1000  gal)  X 0.024 
Design  flow  (MGD) 

Design  surface  settling  rate  (GPD/Ft2) 
If  not  readily  available,  calculate 

SSR  - Design  flow  (MGD)  X 106 
Total  surface  area  (Ft2) 


Columns  49-56:  Design  weir  overflow  rate  (GPD/Ft) 

If  not  readily  available,  calculate 
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WOR  - Design  flow  (MGD)  X 106 
Total  weir  length  (Ft) 

Columns  57-64:  Design  solids  loading  (Lb  SS/Ft2-Day) 

Card  5 - General  Operations 

Columns  1-8:  Peak  flow  (influent)  (MGD) 

Columns  9-16:  Total  daily  flow  (influent)  (MG) 

Columns  17-24:  BOD  (influent)  (mg/1) 

Columns  25-32:  Suspended  solids  (influent)  (mg/1) 

Columns  33-40:  BOD  (effluent)  (mg/1) 

Columns  41-48:  Suspended  solids  (effluent)  (mg/1) 

Card  6 - Process  Operations 

Columns  1-8:  Primary  clarifier  effluent  BOD  (mg/1) 

Columns  9-16:  Primary  effluent  suspended  solids  (mg/1) 

Columns  17-24:  Trickling  filter  effluent  BOD  (mg/1) 
Columns  25-32:  Filter  effluent  suspended  solids  (mg/1) 
Columns  33-40:  Total  recycle  to  filter  (MGD) 

Columns  41-48:  Total  recycle  from  filter  effluent  (MGD) 
Columns  49-56:  Recycle  from  secondary  clarifier  (MGD) 

Activated  Sludge  Type  Plant 
Card  1 

Column  1:  Plant  type;  2 for  activated  sludge  type  plant. 

Columns  2-41:  Same  as  for  trickling  filter  type  plant. 

Card  2 - Primary  Clarifier 

Same  as  for  trickling  filter  type  plant. 

Card  3 - Aeration  Tanks 

Columns  1-8:  Total  tank  volume  (1000  Gal) 

(2  tanks  @ 50,000  gal  ea.  • 100.) 


L 
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Columns  41-48:  Recycle  sludge  volume  (MGD) 

Columns  49-56:  Recycle  sludge  suspended  solids  (mg/1) 

Columns  57-64:  Wasted  sludge  volume  (GPD) 

Columns  65-72:  Average  aeration  basin  disolved  oxygen  (mg/1) 
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IV.  AUTOMATED  INSTRUMENTATION 

A Fisher  Scientific  Accumet  Model  520  Digital  pH/Ion  meter,  a Yellow  Springs 
Model  33  S-C-T  meter,  A Talos  5I4B  graphic  tablet,  and  a Digital  Equipment  Cor- 
poration RT02-BA  data  entry  terminal  were  interfaced  to  the  data  handling  system 
via  a Digital  Equipment  Corporation  PDM-70  programmable  data  mover.  Figure  XI 
is  a block  diagram  of  the  automated  instrumentation  showing  the  inter-connection 
of  these  devices.  Each  of  these  devices  can  operate  in  a stand-alone  mode  inde- 
pendent of  the  data  handling  system  or  under  the  control  of  the  data  handling 
system.  The  Automated  Instrument  Users'  Manual  and  the  Automated  Instrument 
System  Programmer's  Manual  have  been  written  to  provide  a detailed  description 
of  the  operation  of  this  system.  These  Manuals  are  Appendix  V and  VI  of  this 
report.  The  presentation  which  follows  is  intended  to  provide  an  overview  of 
the  automated  instrumentation. 

The  Talos  digitizer  provides  the  capability  to  the  Army  personnel  of 
automating  the  data  reduction  of  field  strip  chart  recordings.  By  placing  a 
strip  chart  on  a )4  inch  by  14  inch  active  area  on  the  Talos  digitizer  surface, 
the  user  through  either  a four  button  cursor  or  a stylus  may  determine  the  x-y 
coordinates  of  any  point  on  this  graph  to  an  accuracy  of  .01  inches.  The  x-y 
coordinates  of  this  point  are  displayed  on  a front  panel  display  and  can  be 
transferred  to  the  data  handling  system  via  the  PDM-70  by  pressing  a cursor 
button  or  touching  the  point  with  the  stylus.  Some  status  information  such  as 
the  cui sor  push  button  is  also  interfaced  to  the  PDM-70.  The  data  handling 
system  can  then  be  programmed  with  the  appropriate  data  reduction  algorithm  to 
arrive  at  the  result.  The  interface  between  the  Talos  digitizer  and  the  PDM-70 
was  fabricated  at  Clemson  University. 

The  Fisher  Scientific  pH/Ion  meter  provides  a 5 digit  display  on  the  front 
panel  of  the  meter  and  is  interfaced  to  the  PDM-70. 
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The  interface  between  the  PDM-70  and  the  pH  meter  was  implemented  at  Clemson 
University  and  contains  a button  with  a light  for  data  transfer.  The  measurement 
value  being  displayed  on  the  front  panel  display  along  with  the  front  panel 
switch  settings  and  over  range  indicator  is  transferred  to  the  PDM-70  by  pressing 
the  lighted  button  on  the  PDM-70  interface. 

The  Yellow  Springs  S-C-T  meter  was  slightly  modified  to  permit  the  meter  to 
give  two  analog  outputs.  One  output  is  proportional  to  the  S-C-T  meter  reading 
while  the  other  output  is  coded  to  give  the  front  panel  switch  settings.  These 
analog  outputs  are  interfaced  to  the  PDM-70  with  a meter  reading  being  trans- 
ferred whenever  a lighted  button  associated  with  the  PDM-70/S-C-T  interface  is 
pressed  by  the  user.  This  interface  was  fabricated  at  Clemson  University. 

Each  measurement  handled  by  the  data  handling  system  must  have  an  ID  and 
a measurement  value.  The  DEC  RT-02-BA  provides  the  user  of  the  automated 
instrumentation  the  ability  to  enter  the  ID  information  via  the  keyboard  for  each 
of  the  samples  whose  value  is  being  measured  by  either  the  interfaced  pH  and 
specific  conductivity  instrumentation.  The  DEC  RT-02-BA  also  has  a 32  character 
alphanumeric  display  and  provides  the  user  with  messages  from  the  data  handling 
system. 

The  PDM-70  programmable  data  mover  provides  a communications  link  and 
formatting  facilities  for  concentrating  and  transferring  data  between  a selected 
source  and  destination.  The  data  may  be  in  an  analog,  parallel  BCD  or  binary, 
and  serial  format.  The  PDM-70  can  be  operated  as  a stand-alone  unit  or  under  the  control 
of  a host  computer.  Integration  into  a computer-based  system  is  easily  accom- 
plished through  a standard  asychronous  serial  interface  on  the  computer. 

The  PDM-70  system  accomodates  both  source  and  destination  modules.  Any 
source  module  under  the  control  of  the  PDM-70  can  communicate  with  any  or  all 
destination  modules.  Source  modules  present  in  the  system  are  as  follows: 
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2 PDM70-D  Eight  digit  BCD  input 

1 PDM70-F  Four  channel  analog  input 

1 PDM70-JR  Serial  I/O  module 

2 PDM70-J  Serial  I/O  module 

The  destination  modules  present  in  the  system  are  as  follows: 

1 PDP70-E  Eight  digit  BCD  output 

1 PDP70-H  Two  channel  analog  output 

2 PDM70-J  Serial  I/O  module 

1 PDM70-JR  Serial  I/O  module 

Communication  or  information  transfer  from  a selected  source  to  the  selected 
destination  modules  is  controlled  by  the  stored  program  being  executed  by  the 
PDM-70.  This  program  is  normally  entered  via  a PDM70-JR  serial  I/O  module  with 
a connected  Teletype  or  remote  computer.  The  programming  commands  are  a sequence 
of  ASCII  characters. 

A user  function  for  OS/8  BASIC  has  been  provided  to  allow  any  BASIC  program 
to  remotely  program  the  PDM-70  and  to  transfer  data  between  the  BASIC  program 
and  the  PDM-70.  In  addition,  the  data  entry  software  for  the  data  handling 
system  supports  the  automated  instrumentation. 

The  OS/8  BASIC  program  WTHDL.BA  is  provided  for  handling  the  Talos  digiti- 
zer data.  This  program  receives  x-y  coordinate  information  from  the  digitizer 
on-line,  converts  these  coordinates  from  the  digitizer  units  to  the  units  of  the 
strip  chart  recording,  and  places  these  converted  coordinates  in  a file.  The 
program  provides  the  user  with  the  capability  to  edit  the  data  received  from  the 
digitizer  and  the  data  within  the  data  files  created  by  this  program.  These  data 


files  can  be  accessed  by  any  OS/8  BASIC  program. 
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V.  DATA  HANDLING 

A data  handling  system  has  been  Implemented  to  automate  the  handling  of 
field  data  during  a water  quality  survey.  Figure  XII  is  a block  diagram  of 
the  data  handling  system  showing  the  Inter-connection  of  the  system  components. 
Table  111  lists  the  subsystem  components  which  comprise  this  system.  Each 
measurement  handled  by  the  system  is  represented  by  6 primary  ID  tags,  5 secon- 
dary ID  tags,  and  a measurement  value  with  each  tag  being  a variable  length 
character  string  Functions  performed  by  the  system  include  data  entry,  report 
generation,  file  management,  analysis,  and  sorts.  The  data  entry  function  is 
capable  of  accepting  one  or  more  raw  measurement  values  for  each  measurement 
and  then  perform  the  appropriate  data  reduction  to  arrive  at  a final  measurement 
value.  The  Data  Handling  System  Users'  Manual  and  System  Programmer's  Manual 
have  been  written  to  provide  a detailed  description  of  the  operation  of  this 
system.  The  presentation  which  follows  is  intended  to  provide  an  overview  of 
the  data  handling  system. 

The  data  handling  system  may  be  located  at  a fixed  remote  site  and/or  may 
be  found  on-site.  A remote  site  would  require  a voice  grade  telephone  line 
data  link  with  the  field  survey  whenever  communications  between  the  data 
handling  equipment  at  the  field  site  and  the  data  handling  system  are  occurring. 
Occassional  assistance  at  the  data  handling  system  site  will  be  required  to 
perform  functions  such  as  to  mount  magnetic  tapes  and  to  start  the  data  handling 
system. 

The  data  handling  system  may  be  viewed  as  a collection  of  programs  which 

I 
I 
I 
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operate  within  the  supervision  of  the  OS/8  executive  on  a DEC  PDP-8/E  minicom- 
puter. The  OS/8  is  a software  system  which  Is  supported  by  a large  variety  of 
hardware  configurations  with  a DEC  PDP-8/E  CPU. 
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A.  OS/8  Minicomputer  System 


89 


PDP-8/E-AE  DEC  8k  Minicomputer  with  Asynchronous  Serial  Interface 

TD8-EM  Dual  DECtape  Transports  with  Control 

MR8-EC  OS/8  ROM  for  TD8-EM 

KL8-EJ  Serial  Interface^for  Mark  Reader 

ASR-33  Teletype 

B.  Hewlett-Packard  7260A  Optical  Mark  Reader 

C.  DEC  PDM-70  Programmable  Data  Mover  System 

PDM70-C  Mounting  Box  with  Power  Supply 
PDM70-P  Programmable  Bus  Control 
PDM70-D  Eight  Digit  BCD  Input  (Digitizer) 

PDM70-D  Eight  Digit  BCD  Input  (pH) 

PDM70-E  Eight  Digit  BCD  Output 

PDM70-F  Four  Channel  Analog  Input  (Conductivity) 

' PDM70-J  Serial  Interface  (Manual  Data  Entry) 

PDM70-J  Serial  Interface  (Optical  Reader) 

PDM70-JR  Serial  Interface  (OS/8  System) 

PDM70-H  Two  Channel  Analog  Output 

D.  Automated  Data  Entry  Equipment 

Fisher  Scientific  Accumet  Model  520  Digital  pH/Ion  Meter  with  Accessories 
Yellow  Springs  Model  33  S-C-T  Meter  with  Accessories 
DEC  RT02-BA  Data  Entry  Terminal 

Talos  Systems  5 1 Graphic  Tablet  with  BCD  and  Four  Positions  Cursor  Options 

E.  DEC  TM8-FA  7"Track  Industry  Compatible  Magnetic  Tape 

F.  Other 

3 - Omnitec  701-B  Originate  Couplers 
DEC  KL8-M  Modern  Control  Interface 
DEC  BE8-A  Omnibus  Expansion  (TM8-FA) 


TABLE  III  DATA  HANDLING  SYSTEM  SUBSYSTEM  COMPONENTS 
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Each  of  these  programs  operate  in  conjunction  with  data  files  found  on  the  mini- 
computer's mass  storage  devices.  Because  of  core  requirements,  a particular 
data  handling  operation  is  often  implemented  by  chaining  several  of  these  pro- 
grams together.  Specfal  user  programs  can  be  readily  implemented  using  languages 
such  as  BASIC,  FORTRAN  II,  and  FORTRAN  IV,  to  analyze  the  information  found 
within  the  data  base  created  by  the  data  handling  system  software. 

All  measurements  are  represented  by  6 primary  ID  tags,  5 secondary  ID 
tags,  and  a measurement  value  with  each  ID  tag  being  a variable  length  charac- 
ter string.  The  primary  ID  tags  uniquely  identify  which  measurement  while  the 
secondary  ID  tags  give  information  about  the  measurement.  Figures  XIII,  XIV, 
and  XV  respectively  show  the  measurement  format,  primary  ID  tag  assignments, 
and  secondary  ID  tag  assignments.  More  than  one  quality  ID  may  be  used  to 
describe  each  measurement  with  commas  separating  these  tags.  The  first  example 
measurement  in  Figure  XVI  indicates  that  conductivity  was  measured  on  day  1 
of  the  survey  at  sample  point  Cl  with  this  measurement  being  the  second  discrete 
measuring  during  period  2.  The  water  sample  was  turbid,  light  color,  and  hot 
in  appearance.  The  analysis  was  performed  by  Drake  using  method  1 to  yield 
a value  of  700  and  the  analysis  yielded  an  accurate  result. 

Data  and  commands  are  usually  entered  into  the  system  on-line  via  either 
the  Teletype  keyboard,  a Hewlett-Packard  optical  mark  reader,  the  Talos  graphic 
digitizer,  the  DEC  manual  data  entry  station,  Yellow  Springs  S-C-T  meter,  or 
Fisher  pH/Ion  meter.  The  digitizer,  manual  data  entry  station,  specific 
conductivity  Instrumentation,  and  pH  instrumentation  are  interfaced  to  the 
system  via  the  DEC  PDM-70  programmable  data  mover.  Then  data  when  received 
is  extensively  checked  to  Insure  that  only  valid  data  Is  entered  into  the  data 
files  with  numerous  error  messages  being  generated  to  Identify  data  entry 
errors.  The  system  is  capable  of  accepting  one  or  more  raw  measurement  values 
for  each  measurement  on  Input  and  then  perform  the  appropriate  data  reduction 
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to  arrive  at  a final  measurement  value.  The  particular  data  reduction  algorithm 
selected  Is  determined  by  the  ID  tags  associated  with  the  measurement.  Commands 
are  provided  for  the  user  to  enter  constants  relative  to  each  data  reduction  al- 
gorithm and  to  select  one  of  several  Input  options.  When  the  optical  mark  reader 
is  used  on  Input,  a special  mark  card  Is  designed  for  each  application.  This 
card  has  fields  for  command  Information,  ID  Information,  and  value  entry.  The  use 
of  mark  cards  permits  users  of  the  system  to  enter  data  on  a card  with  any  combina- 
tion of  pencil  marks  and  punches.  Mark  sense  cards  provide  an  effective  low  cost 
data  entry  method  especially  when  this  data  Is  belhg  entered  on-line  by  several 

system  users.  A unique  translate  table  Is  provided  for  each  column  to  translate 
each  character  received  from  this  reader  Into  the  appropriate  character  string. 

The  resulting  character  string  which  results  from  each  card  or  sequence  of  cards 
becomes  the  measurement  input  accepted  by  the  system. 

The  file  management  software  gives  the  user  of  the  system  the  capability  of 
adding,  deleting,  or  changing  measurements  within  a specified  data  file.  In  addi- 
tion the  system  can  transfer  files  between  devices,  merge  and  delete  files,  and 
list,  zero,  and  compress  directories  on  the  mass  storage  devices. 

An  important  system  function  is  the  capability  to  generate  reports  of  selected 
data  on  a designated  output  device  using  a wide  variety  of  formats.  Report  genera- 
tion is  accomplished  in  three  Job  steps  which  are  report  specification,  sort,  and 
report  preparation.  The  general  report  format  Is  determined  by  the  order  in  which 
the  measurements  have  been  sorted  and  ts  shown  in  Figure  XVII.  Note  that  each 
page  of  the  report  shows  a tabulation  of  measurement  values  with  three  primary  ID 
tags  being  held  constant.  The  comment  ID  tag,  ID6,  when  nonzero  qualifies  each 


measurement  value  with  a symbol  to  Indicate  that  the  actual  measurement  is  less 
than,  greater  than,  or  approximately  equal  to  the  value  shown. 
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f D!  - Day  ID 
ID2  - Point  ID 
IDi  - Period  ID 
IDA  ■ Parameter  ID 
ID5  ■ Discrete  ID 
ID6  * Comment  ID 

FIGURE  XIV  PRIMARY  ID  TAG  ASSIGNMENTS 


1 07  - Quality  ID 
Turbid  (T) 

Clear  (CL) 

Light  Color  (LC) 

Intense  Color  (1C) 

Hot  (H) 

Cold  (C) 

Greasy  (G) 

Nonhomogeneous  (NH) 
Imperfect  Preservation  (IP) 
Questional  Stability  (QS) 

ID8  - Accuracy  ID 

Accurate  Result  (AR) 

10*  (10*) 

Order  of  Magnitude  (0M) 
Inaccurate  Results  (IR) 

ID9  - Chemist  ID 
I DIO  - Method  ID 
I D1 1 • Unused 


FI  BURE  XV  SECONDARY  ID  TAG  ASSI6NMENTS 


1 Cl  2 COND  2 


T,LC,H  AR  DRAKE  Ml  700 
DRAKE  7.4 


10  SI  3 PH 
2 15  1 TOC  25  BROWN  200 
FIGURE  XVI  EXAMPLES  OF  MEASUREMENTS 


Second  Order 
Fourth  Order 


Heading 
First  Order 

Third  Order 

Fifth  Order 


Footnotes 


FIGURE  XVII  REPORT  FORMAT 


94 


The  report  spec  I f icat Ion  step,  implemented  with  OS/8  BASIC,  creates  the 
file  INFO: INFO. IN  based  on  Information  received  from  the  operator  of  the  sys- 
tem. This  file  contained  one  record  for  each  report  requested  during  the  re- 
port specification  step.  Each  record  contains  Information  such  as  report 
format  (order  of  sort),  options  (report  or  print-plot),  ID  tag  lists,  and 
segment  size  (number  of  columns).  The  order  in  which  the  ID  tags  appear  with- 
in the  lists  determine  the  order  they  appear  in  the  report. 

At  the  completion  of  the  report  generation  step,  the  sort  software, 
implemented  with  PAl8,  is  automatically  loaded  by  the  system.  The  specified 
data  file  Is  searched  sequentially  to  locate  and  retrieve  those  measurements 
which  satisfy  the  ID  lists  within  the  file  INFO. IN.  These  measurements  are 
placed  within  a core  buffer.  The  measurements  within  this  buffer  are  sorted 
according  to  a specified  order,  order  of  the  ID  tags  within  the  ID  list,  and 
segmentation  size  and  then  places  these  measurements  within  a designated  tem- 
porary file.  A temporary  file  is  created  for  each  report.  Measurements  with 
duplicate  ID  tags  are  detected  and  placed  within  a duplicate  data.fi le. 

The  sort  step  generally  utilizes  only  the  primary  ID  tags  while  sorting 
the  measurements.  The  secondary  ID  tage  are  either  ignored  or  result  in 
footnotes  on  the  reports  generated. 

The  report  preparation  step,  implemented  with  OS/8  BASIC,  utilizes  the 
file  INFO. IN,  the  system  definition  file  SYSDEF.AR,  and  the  temporary  data  files 
to  produce  on  the  output  device  the  designated  reports.  The  system  definition 
file  is  used  to  convert  the  Internal  ID  tag  formats  to  an  external  format  re- 
quired by  the  report.  As  each  ID  tag  has  several  external  representations, 
the  particular  representation  depends  upon  the  usage  within  the  report.  For 
example,  a 24  character  external  representation  Is  used  to  label  a row  of  the 
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report  while  a 6 character  representation  labels  columns.  The  system  defini- 
tion file  specifies  the  number  of  digits  following  the  decimal  point  to  use 
when  printing  each  numerical  value  based  on  the  ID<*  tag.  Also  found  are  spe- 
cifications regarding  margins,  footnotes,  and  header  Information. 

The  report  preparation  step  was  specifically  designed  to  prepare  each  page 
of  the  report  in  a format  which  can  be  reproduced  directly  In  a final  report. 
Each  page  is  generally  printed  on  a 8 1/2  x 11  page  with  standard  margins. 

The  system  has  a DEC  TM8-FA  7*track  Industry  compatible  magnetic  tape 
unit  to  provide  data  transfer  capability  between  this  system  and  the  Army 
UNIVAC  1108.  The  subroutine  NTRAN (FORTRAN  callable)  provides  FORTRAN  V pro- 
grams written  for  the  Army  UNIVAC  1108  with  the  capability  for  reading  and 
writing  magnetic  tapes  compatible  with  the  OS/8  system.  The  subroutine  NTRAN 
is  well  documented  In  the  UNIVAC  FORTRAN  V Library  Programmer  Reference  Manual 
(UP-7876).  The  OS/8  magnetic  tape  format,  required  by  NTRAN,  Is  found  in  the 
System  Programmer's  Manual  for  the  data  handling  system. 

The  performance  of  the  system  Is  now  highly  dependent  on  the  10  character- 
per-second  (cps)  printer  rate  of  the  Teletype.  This  data  rate  requires  the 
system  to  be  committed  for  large  blocks  of  time  to  the  printing  of  reports, 
messages,  and  plots.  The  replacement  of  the  Teletype  with  a Digital  Equipment 
Corporation  LA36  terminal  Is  recommended  and  would  Increase  the  printing  rate 
to  30  cps  and  provide  a 132  print  positions  per  line.  The  rental  of  a 30  cps 
terminal  during  the  third  year  from  Western  Union  Data  Services  Company  demon- 
strated the  performance  Improvements  which  can  be  obtained. 

The  addition  of  the  second  hard  copy  device  such  as  the  DEC  LA36  terminal 
would  permit  one  hard  copy  device  at  the  field  site  and  the  other  hard  copy 
device  to  be  at  the  OS/8  site.  The  use  of  two  hard  copy  devices  would  also 


permit  operation  of  the  OS/8  system  from  either  location. 

Several  system  tables  are  used  by  the  system  software  to  define  a parti- 
cular data  handling  application  to  the  system.  The  system  definition  file 
(SYSDEF.AR)  provides  Information  which  includes  lists  of  valid  ID  tags,  upper 
and  lower  limits  for  measurement  values,  data  fomats,  and  various  heading  and 
column  information  used  during  report  generation.  The  translation  table  file 
(TRHTBL.AR)  contains  the  translation  tables  for  the  mark  sense  cards.  The 
method  file  (METHOD. AR)  specifies  for  each  parameter  the  particular  data  con- 
version algorithm  within  the  function  file  (FNCTN.AR)  to  use  on  input. 

Version  1 1 1 of  the  data  handling  system  was  successfully  demonstrated  by 
handling  survey  data  generated  by  the  USAEHA  during  the  Letterkenny  Army 
Depot  survey  in  August  1971*.  This  version  of  the  system  did  not  include  any 
of  the  third  year  improvements.  The  system  accepted  raw  measurement  data  for 
each  parameter  from  the  chemist  and  performed  the  necessary  data  reduction 
to  arrive  at  a final  measurement  value.  The  various  calibration  curves  used 
by  the  system  were  determined  from  calibration  data  entered  into  the  system 
on  a dally  basis.  Data  validation  was  accomplished  by  having  the  system  pro- 
vide each  chemist  during  data  entry  with  a report  showing  calibration  data, 
raw  measurement  data,  and  the  computed  final  measurement  result.  The  system 
was  used  an  average  of  approximately  two  hours  each  day  for  all  data  handling 
functions  during  the  survey. 

Version  IV  of  the  data  handling  system  satisfies  the  third  year  contract 
objectives  and  differs  from  the  version  1 1 1 of  the  system.  Version  IV  in- 
cludes the  following  Improvements  not  found  in  version  III: 

1.  new  measurement  format  (secondary  ID  tags) 

2.  new  mark  card  format 


3. 


improved  data  reduction  for  the  chemists 
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4. 

Interfaced  pH  and  specific  conductivity  Instrumentation 

5. 

improved  software 

6. 

industry  compatible  magnetic 

tape 

7. 

8K  core  expansion  to  16K 

8. 

interfaced  graphic  digitizer 

9. 

interfaced  manual  data  entry 

station 

This  version  of  the  system  satisfies  the  data  handling  objectives  established 
on  September  12,  1974  in  a conference  between  the  Army  and  Clemson  University 
personnel.  This  version  of  the  system  has  not  been  evaluated  in  a field  en- 
v i ronment . 

The  two  DEC  tape  units  provide  all  the  file  storage  for  the  data  handling 
software.  As  most  of  the  software  both  access  and  update  the  data  base  found 
within  several  files,  the  execution  speed  of  this  software  is  highly  dependent 
on  the  average  access  time  of  these  two  DEC  tape  units.  The  addition  of  a 
disk  to  the  system  would  significantly  Increase  the  execution  speed  of  the 
software  by  substantially  reducing  the  time  required  to  access  the  disk  data 
files.  However,  the  speed  of  the  current  system  configuration  adequately 
handled  the  Letterkenny  Depot  data  handling  demonstration. 


